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^ ABSTRACT 

We present VLT/VIMOS-IFU emission-line spectroscopy of a volume limited sample 
of 18 southern ULIRGs selected with z < 0.09 and S < 10. By covering a wide range of 
ULIRG types, including many systems that have received very little previous attention, 
this dataset provides an important set of templates for comparison with high-redshift 
Q^J galaxies. We employed an automated Gaussian line fitting program to decompose the 

emission line profiles of Ha, [N il], [S n], and [O i] into individual components, and chart 
the Ha kinematics, and the ionized gas excitations and densities. 11/18 of our galaxies 
show evidence for outflowing warm ionized gas with speeds between 500 and a few 
1000 kms -1 , with the fastest outflows associated with systems that contain an AGN. 
Our spatially resolved spectroscopy has allowed us to map the outflows, and in some 
cases determine for the first time to which nucleus the wind is associated. In three of 
our targets we find line components with widths >2000 kms -1 over spatially extended 
regions in both the recombination and forbidden lines; in two of these three, they are 
^— H associated with a known Sy2 nucleus. Eight galaxies have clear rotating gaseous disks, 

^ and for these we measure rotation velocities, virial masses, and calculate Toomre Q 

CO parameters. We find radial gradients in the emission line ratios in a significant number 

of systems in our study. We attribute these gradients to changes in ionizing radiation 
field strength, most likely due to an increasing contribution of shocks with radius. 
We conclude with a detailed discussion of the results for each individual system, with 
reference to the existing literature. 

Our observations demonstrate that the complexity of the kinematics and gas prop- 
C^l erties in ULIRGs can only be disentangled with high sensitivity, spatially resolved IFU 

observations. Many of our targets are ideal candidates for future high spatial resolution 
follow-up observations. 

k> Key words: galaxies: interactions - galaxies: ISM - galaxies: starburst - galaxies: 

active - ISM: kinematics and dynamics - ISM: jets and outflows. 
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1 INTRODUCTION 

Ultra-Luminous InfraRed Galaxies (ULIRGs) are the most 
IR luminous population in the local universe (Sanders & 
Mirabel 1996). While rare locally, the appearance of large 
numbers of luminous objects in far-IR/submm blank field 
surveys (e.g |Smail et al.|[l 997i Hu ghes et al.||1998| |Barger| 
et al.||1998| |Eales et al.||2000| |Coppin et al.||2006| ) suggests 



* Based on observations collected at the European Organisation 
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der programmes 080.B-0085 and 383.B-0372. 
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that they were more important in the high-z universe. The 
process that triggers ULIRG activity appears to be a major 
merger between two disk gala xies (|Clements et al.|1996| |Duc| 



et al.|1997| |Borne et al.|1999[ |2000) , which consequently ini- 
tiates a (at least partially heavily obscured) starburst and/or 
AGN. In fact ULIRGs are implicated as hosting ^50% the 
star formation in the universe at z ~ 1 a nd above (|Perez- 



Gonzalez et al.|2005| |Le Floc'h etal. 2005 ; B ethermin et al. 



2011 Clements et al. in prep.). Since most massive galaxies 
are thought to have experienced at least one major merger 
(Ma n et al.|2011 ), it is therefore essential to study the phys- 
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ical mechanisms at work in local ULIRGs where the best 
spatial resolution is achievable. 

Hydrodynamic simulations predict that as galaxies col- 
lide and coalesce, gravitational forces and dissipative shocks 
funnel gas toward the centre, triggering a nuclear starburst 



and fuelling the growth of the central black hole ( Barnes & 
Hernquistfl992| |Cox et al.||2008| ). The end product of such 
a "wet" (i.e. gas-rich) merger is likely to be an elliptical 



galaxy, as predicted by theory (Barnes &; Hernquist 1996) 



and found by observations of local ULIRGs (e.g., Wright 



et al. 1990 Genzel et al. 2001 Dasyra et al. 2006b ). The high 



IR luminosities of ULIRGs therefore arise through the dust 
re-processing of UV radiation from the concomitant star- 
burst/AGN. Recent studies, however, are challenging this 
idea 



2010 



(e.g. : 



Forster Schreiber et al. 2009 Chien & Barnes 



Cisternas et al.||2011 ), and the emerging view is that 



merger- induced starbursts will often have both a nuclear 
and extended component in the form of massive star- forming 
clumps in the disk ( |Bournaud||2011| . 

Although some ULIRGs are known to host an AGN, 
their role, and the relationship between ULIRGs and 
quasars, is still unclear ( Alexander et aL||2005 ). The preva- 
lence of optically visible AGNs increases with increasing in- 



frared luminosity (Veilleux et al. 1999) and merger stage 



( Veilleux & Rupke 2002 ) , supporting the canonical view that 



the dust obscuring an AGN is removed through feedback 
processes, revealing an optical quasar ( |Sanders et al.|1 988). 
However, this evolutionary sequence has also been brought 
into question by recent work that, for example, shows that 
many ULIRGs are evolving into ellipticals of lower mass 
than known quasar hosts (|Genzel et al.||200l) |Veilleux fe| 
Rupke| [20021 |Tacconi et al.||2002| |Dasyra et al.| |20Q6a|b[ 
2007). Evidence for feedback that is clearly AGN-driven in 
LIRGs or ULIRGs has thus far re mained elusive ([Ru pke 



eTaLl[200 5a; Rup ke fe Veilleux||2011| |Muller-Sanchez et al. 
20TT| 



Nevertheless, winds with significant mass and energy 



flux have been observed in ULIRGs (Clements et al.|2002 



Martini [20051 [2006] |Rupke et alTl|2005b|a| |Fu & Stockton| 
2009). Theory predicts that this feedback regulates the 
growth of supermassive black holes and their accompany- 



ing starbursts (Di Matteo et al.||2005 | Veilleux et al.||2005 



Narayanan et al.|2 008; Hop kins &; Elvis|2010 ), making them 
an indispensable element in models of galaxy evolution. 
However, exactly how this feedback regulation occurs is still 
very poorly understood. 

Warm ionized and cool neutral gas is thought to re- 
side in clouds advected into the outflow at the shear in- 
terface between the hot wind and the quiescent disk gas 



(e.g. Tenorio-Tagle Sz Muhoz-Tunon 1998). Emission line 



observations have the ability to trace this warm ionized 
phase, and studies of local starburst ing galaxies have shown 
that superwinds are common (e.g., |Heckman et al.||1990| 



(mostly long-slit) studies of the cool gas component via the 
Nai AA5890,5896 interstellar absorption line. These studies 
have found that ^80% of ULIRG galaxies show evidence for 
neutral outflows with speeds of 300-700 kins" 1 and high 
mass outflow rates of order 10-100 M yr _1 that appear 
to scale with Lir, SFR and galaxy mass (|Heckman et al.| 
'20001 |Rupke et al.|2002| |Martin|2005l |Rupke et al.|2005b|a| 



Martin 2006) 



The major limitation of such absorption line studies is 
the limited radial extent over which the absorption lines can 
be probed. At large radii the surface brightness of the con- 
tinuum source (starlight) becomes too faint to probe any 
absorbing material, whereas ionized gas emission is inde- 
pendent of the stellar surface brightness and can often be 
observed to much larger distances. Thus the last decade has 
seen a steady flux of spatially resolved optical emission line 
studies of (U) LIRGs, both on a detailed case-by-case basis 
(e. g.,|Colina et al.|2000]|Murphy et al.|2001||Arribas fe CoT 
ina|2003| ILipari et al.||2004a|b| |Bedregal et al.|2009| |Bendo 



et al.| |2009i Sh ih fe Rupke||2010| |Rupke fc yeilleux||2011 



and as part of small surveys (Colina et al.| 2005 |Arribas 



et al.) [20081 JGarcia-Marm et al.||2009a|b| |Alonso-Herrero 



et al.||2009[|2010| |Monreal-Ibero et al.||2010| ). To date no 
integral field survey of the emission-line characteristics of 
local ULIRGs has been carried out with spatial resolutions 
better than l" or spectral resolutions better than R^IOOO. 
For these reasons, and to capitalise on the recent interest, 
we have obtained VLT/VIMOS integral field spectroscopy 
of a volume limited sample of 18 southern ULIRGs, with 
spatial and spectral resolutions of ~0f.'7 and R=3100. Our 
sample covers a range of ULIRG types, including starburst 
dominated, LINER, Sy2, starburst +Sy 2, LINER+Sy2, (but 
not the rare Syl class of which there are no \ow-z southern 
hemisphere examples), distant pairs, close pairs, and fully- 
coalesced systems, and has allowed us to measure detailed 
gas dynamics, excitations and densities in a large sample of 
local ULIRG systems for the first time. This paper forms 
part of an ongoing programme to study the gas dynamics in 
ULIRGs; a detailed investigation of IRAS 19254-7245 ("The 
Superantennae" ) was the subject of our first paper (Bendo 



|et al.||2009[ hereafter [Paper l\ . 

The remainder of this article is organised as follows: In 
Section [2] we briefly describe the sample selection, data ac- 
quisition and reduction, and describe how we modelled the 
emission line shapes using multi-Gaussian fitting. In Sec- 
tion [3] we present our emission line maps together with line 
diagnostic plots, and describe a number of measured and 
derived properties of the galaxies in the whole sample (Sec- 
tion |3.1| ) . We then discuss our results with respect to the 
gas dynamics, line ratios and evidence for outflows in Sec- 
tion [4] and summarise our main conclusions in Section [5] In 
the Appendix we give a detailed case-by-case description of 
the emission line maps for each target, with reference to the 
existing literature. 



Marl owe et"aLl[l99 5 Martin 1998 Meurer 2004). Detailed 
studies of individual nearby cases highlight their complex- 
ity, both morphologically and kinematically (e.g., Westmo- 
quette et al.|2007[ [2009 2011 ). Absorption- line spectroscopy, 
however, has proved a much easier method for identify- 
ing and studying outflows since this technique relies on the 
background continuum light from the galaxy, which is much 
brighter than the emission lines. Up to now much of our 
knowledge about ULIRG outflows has therefore come from vey of ULIRGs (Kim &; Sanders 1998). Observations were 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Data acquisition and reduction 

Our sample was selected from the the IRAS 1 Jy sur- 
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performed with the Visible MultiObject Spectrograph (VI- 
MOS; [LeFevre et al.|[2003| ) integral field unit (IFU) at the 
Very Large Telescope in 2007 (080.B-0085) and 2009 (383.B- 
0372). A detailed description of the data reduction proce- 
dure was given in |Paper I[ so we only give a brief recap 
here. The data were selected to be nearby (z < 0.09) and 
accessible from the VLT (5 < 10). They were taken using 
the HR red erism (6300-8700 A) giving a spectral resolution 
of R=3100, and with the / . / 67 spaxel -1 spatial resolution. 
A 5-point dither pattern was used, which provided redun- 
dancy within the central region and ensured the availability 
of blank sky for background measurements. Four exposures 
were taken at each dither position; the total exposure times 
are given in Table ^ The seeing during the observations was 
~0.7-l arcsec. At the time of observation, two of the com- 
plete sample of ULIRGs, IRAS 05189-2524 and IRAS 06035- 
7102, had already been observed with the VIMOS IFU as 



part of the |Arribas et al. ( 2008 ) study. The reduced dat- 
acubes for these galaxies were kindly made available to us 
for inclusion in this study (S. Arribas, private communica- 
tionfl 

Basic data reduction was performed following the 
standard ESO reduction pipeline tasks implemented in 
GASGANO. For each pointing and for each quadrant, the me- 
dian background spectrum was calculated and subtracted 
from the data. Following this, the spectra from each point- 
ing were mapped into individual spectral cubes, and then 
the cubes were median combined to produce the final spec- 
tral cube. An example full HR Red spectrum from IRAS 
00335-2732 is shown in Fig. [I] with the emission lines exam- 
ined in this paper labelled. Finally, a S/N cut was applied to 
each cube based on the integrated line flux in the Ha- [Nil] 
region to remove spaxels with no detected line emission. 



2.2 Decomposing the emission line profile shapes 



Following the methodology first presented by Westmoq uette| 
et al. ( |2007| , we fitted multiple Gaussian profile models to 
the Ha, [Nil], [S n] and [O i] lines using an IDL \ 2 fitting 
package called pan (Peak ANalysis; Dimeo 2005; West mo- 
quette et al.|2007| to quantify the gas properties observed in 
each IFU field. The high S/N and spectral resolution of these 
data have allowed us to quantify the line profile shapes of 
these lines to a high degree of accuracy. Our fitting method- 
ology differs slightly to that employed in |Paper I| since we 
have had to optimise the process for the automatic multi- 
Gaussian decomposition of the emission lines in our entire 
sample. 

Each line in each spectrum was fitted using a single, 
double, and triple Gaussian component initial guess (no 
lines were found to need more than three components) . Line 
fluxes were constrained to be positive and widths to be 
greater than the instrumental contribution to guard against 
spurious results. We made no attempt to fit an Ha stellar 
absorption component since it was never needed to obtain a 
good fit. 

The Ha and [Nn]A6583 lines were fit simultaneously, 
where for each component the wavelength difference between 



1 Note the observations in our programme are a factor of four 
deeper than those in the study of Arribas et al. ( 2008| >. 



the two Gaussian models was constrained to be equal to 
the laboratory difference, and the FWHMs equal to one an- 
other. This constrained approach dramatically improves the 
quality of the fits, particularly for spectra with lower S/N. 
Separately, we fitted the [Sn] doublet also simultaneously, 
and the [OI] line individually. Multi-component fits were 
run several times with different initial guess configurations 
(widths and wavelengths) in order to account for the varied 
profile shapes, and the one with the lowest % 2 fit statistic 
was kept. However, we note that the % 2 minimisation rou- 
tine employed by PAN is very robust with respect to the 
initial guess parameters. 

To determine how many Gaussian components best fit 
an observed profile (one, two or three), we used a likelihood 
ratio test to determine whether a fit of n+1 components 
was more appropriate than an n-component fit. This test 
says that if the ratio of the \ 2 statistics for the two fits falls 
above a certain threshold, then the fits are considered sta- 
tistically distinguishable, and the one with the lower % 2 can 
be selected. Here we determine the threshold ratio by vi- 
sual inspection of a range of spectra and fits. We are aware 
that iProtassov et al. ( 20021) show that using any kind of 



likelihood ratio test (including the F-test) in this kind of 
situation is statistically incorrect, but given the absence of 
a statistically correct alternative that could be applied sen- 
sibly to the volume of data presented here, we have chosen 
to opt for this generalised x 2 ratio test approach. 

This test, however, only tells us which of the fits (sin- 
gle, double or triple component) is most appropriate for the 
corresponding line profile. Experience has taught us that we 
need to apply a number of additional tests to filter out well- 
fit but physically improbable results. We specified that the 
measured FWHM had to be greater than the associated er- 
ror on the FWHM result (a common symptom of a bad fit), 
the fluxes of all components should be >0, and we rejected 
any fits where Xsingie/Xdoubie = (another symptom of a 
bad fit). 

To estimate the uncertainties on the fit parameters, we 
used a bootstrap Monte Carlo technique, assuming that the 
noise is uncorrelated. Using the fit results determined from 
the above step as the initial guess, we re- fitted each profile 
100 times, each time seeding the spectrum with randomly 
generated noise based on the standard deviation of the con- 
tinuum. The standard deviation on each parameter (flux, 
FWHM and radial velocity) over the 100 results was then 
adopted as the associated uncertainty for that parameter, 
while the mean of those results was adopted as the actual 
value. Uncertainties range from 1-30% depending on the 
S/N of the spectrum and how many components are fit. 

The final step of the fitting procedure was to assign 
each fit component to a particular map in such a way as to 
limit the confusion that might arise during analysis of the 
results, such as discontinuous spatial regions arising from 
incorrect component assignments. For triple fits, we specified 
CI to be the brightest and C3 the faintest (regardless of the 
relative velocities of the components), except for the galaxies 
that exhibit a very broad Ha component. For IRAS 05189- 
2524 we specified CI to be the reddest and C3 the bluest 
(regardless of their relative flux), and for IRAS 06035-7102 
we specified C3 to be the broadest then CI the brightest of 
the remaining two components (regardless of their relative 
velocities). And for IRAS 19254-7245, we set C3 to be the 
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Table 1. VIMOS ULIRG galaxy sample. 



Galaxy 


a (J2000), 5 (J2000) 


z a 


D a 

(Mpc) 


\og(L m /L Q ) b 


Sp type c 


SFR d 
(Moyr- 1 ) 


Exposure 
time (s) 


IRAS 00198-7926 


00 21 52.9, 


-79 10 08.0 


0.0728 


300 


12.01 


Sy2 (AGN in N N e ) 


53 


12000 


IRAS 00335-2732 


00 36 00.5, 


-27 15 34.5 


0.0693 


282 


11.93 


SB 


117 


12000 


IRAS 03068-5346 


03 08 21.3, 


-53 35 12.0 


0.0745 


304 


11.96 


SB 


126 


12000 


IRAS 05189-2524 


05 21 50.5, 


-25 21 46.3 


0.0426 


173 


12.19 


Sy2 


80 


3000 


IRAS 06035-7102 


06 05 11.2, 


-71 03 05.5 


0.0795 


325 


12.26 


SB (AGN in N w ? e ) 


182 


2880 


IRAS 09111-1007W 


09 13 38.8, 


-10 19 20.3 


0.0541 


220 


11.98 


SB+Sy2 


66 


12000 


IRAS 12112+0305 


12 13 46.0, 


+02 48 38.0 


0.0733 


300 


12.28 


LINER 


263 


12000 


IRAS 14348-1447 


14 37 38.4, 


-15 00 22.8 


0.0827 


345 


12.30 


LINER (N S w = Sy2-f) 


275 


9600 


IRAS 14378-3651 


14 40 58.9, 


-37 04 33.0 


0.0676 


280 


12.24 


LINER+Sy2 


120 


12000 


IRAS 17208-0014 


17 23 21.9, 


-00 17 00.4 


0.0428 


175 


12.25 


SB (LINERS) 


245 


12000 


IRAS 19254-7245 


19 31 21.6, 


-72 39 21.7 


0.0617 


250 


11.91 


Sy2 


56 


12000 


IRAS 19297-0406 


19 32 20.7, 


-04 00 06.0 


0.0857 


353 


12.36 


SB 


316 


12000 


IRAS 20046-0623 


20 07 19.4, 


-06 14 26.0 


0.0844 


348 


11.97 


SB 


129 


12000 


IRAS 20414-1651 


20 44 18.2, 


-44 33 37.7 


0.0871 


358 


12.26 


SB 


251 


12000 


IRAS 20551-4250 


20 58 27.4, 


-42 39 03.1 


0.0428 


175 


11.98 


Sy2 h 


132 


6000 


IRAS 21504-0628 


21 53 05.5, 


-06 14 49.9 


0.0776 


320 


11.92 


Sy2? 


57 


12000 


IRAS 22491-1808 


22 51 49.3, 


-17 52 23.4 


0.0778 


320 


12.09 


SB 


170 


12000 


IRAS 23128-5919 


23 15 47.0, 


-59 03 16.9 


0.0446 


182 


11.96 


SB+Sy2 


63 


6000 



a Reference: NED. 

b 8-1000 urn infrared fluxes from the IRAS Faint Source Catalog (FSC) version 2 ( Moshir fe et al.||l990| , recalculated with updated 
redshifts where required. 

c References: |Farrah et al.| ( |2003| >, NED, except where otherwise indicated. 

d Star formation rate derived from Ljr and scaled according to AGN input (see text for scaling factors) 
e N n = northern nucleus. Reference: this stud y 
(|l988t;|Nakajima et al.||l99l] ) 



Sanders et al. 



Arribas &; Colina ( 2003 ) re-classify this galaxy as a LINER due to the effect of differential extinction on the optical nuclear spectrum. 



Franceschini et al.| fl2003| ); |Ptak et ah] ( |2003| ); 



Far rah et al. 



(2007) 



broadest component then CI as the narrowest (regardless of 
their relative flux) , and for a few double fits we manually set 
a few very broad components to C3 to make the maps more 
consistent. In all other double-component fits we specified 
CI to be the brightest of the two. 



Many galaxies show 2-3 line components in nuclear re- 
gions. However, in some cases the difference between the 
components in the FWHM and/or radial velocity plots can 
be difficult to see due to the small differences in values com- 
pared to the range of the colour bar. This is an inherent 
limitation in the dynamic range that can be represented in 
these type of plots. To highlight this we have chosen an ex- 
ample from the IRAS 20046-0623 map. In Fig. [2] we show 
the Ha+[N n] line profiles and best-fitting single and double- 
Gaussian models for the spaxel indic ated with a white cross 
in the Ha C2 FWHM map in Fig. 



3m 



Both the x an d 
residuals demonstrate that a two-component fit is the most 
appropriate choice, although on the maps it is difficult to 
see the difference between the two components. 



The emission line maps are shown in Figs. |3a{{3r| Each 
figure includes example Ha+[Nn] line profiles and best- 
fitting Gaussian models from selected spaxel(s) within the 
IFU field, together with the integrated spectrum of the 
whole system. These examples have been chosen to demon- 
strate the quality of the spectra and the accuracy of the 
line- fitting. 



3 RESULTS 



Emission line maps created from the line profile fits de- 
scribed above are shown in Figs. [3aJ|3r] For each system 
we show the HST F814W image (or DSS image if HST 
imaging is not available), and for each identified line com- 
ponent, the Ha flux, FWHM and radial velocity maps, the 
[Nn]A6583/Ha, [S n] (A6717+A6731)/Ha and [Oi]A6300/Ha 
flux ratio maps, and the electron density map derived from 
the [S n] doublet ratio (where we have summed all line com- 
ponents together) using the IRAF nebular task and assum- 
ing T e = 10 4 K. 



[Nn]A6583/Ha vs. [S n] (A6717+A6731)/Ha flux ratios 
are also plotted for each galaxy on a spaxel-to-spaxel basis in 
Figs. |4af[4c| together with predictions from the photoioniza- 
tion and shock models of Dopita et al.| ( [2006| and | Allen et al.| 
(2008), and (where appropriate) the dusty AGN models of 
| Groves et al. (2004). Without the [Olll] line we cannot plot 
any of the traditional BPT diagrams ( |Baldwin et al.|1981 ), 
however the [Nil] /Ha vs. [Sll]/Ha plot can still be instruc- 
tive to investigate the contribution of harder radiation fields 
(AGN/shocks) to the ionization. 



We discuss these maps in the context of the whole sam- 
ple in Section^] and describe them in detail for each galaxy 
individually in Appendix [A] 
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3.1 Measured and derived properties of the 
galaxies in the whole sample 

Integrated SFRs for each of our galaxies are listed in Table[l] 
and are derived from a modified version of the Kennicutt 



(1998) calibration of IR luminosity to SFR: 



SFR: 



5.8 x 10 9 L 







(i) 



where a is a factor equal to the fraction of the IR luminosity 
powered by star formation. Following|Rupke et al. (2005b|a]), 



we assume a — 0.8 for SB ULIRGs, a — 0.3 for Sy2 galaxies 
that show either broad lines in the near-IR or faint broad 
lines in the optical (e.g. IRAS 05189-2524), and a = 0.4 for 
all other Sy2s. 



Following Monreal-Ibero et al. (2010), we morphologi- 



cally classify the galaxies into three classes (based loosely 
on the |Veilleux et al.|2002| scheme), using the HST F814W 
imaging (or DSS imaging where HST imaging is not avail- 
able). Our classifications are shown in Table[2] Class repre- 
sents objects that appear to be single isolated galaxies with 
a relatively symmetric morphology and no obvious evidence 
for strong past or ongoing interaction. Class 1 includes ob- 
jects in a pre-coalescence phase with two nuclei separated 
by a projected distance of D > 1.5 kpc. This limit was 
chosen because theoretical models predict a fast coalescence 



phase after the nuclei become closer than 1.5 kpc (Mihos & 
Hernquist 1996 N aab et aL||2006 ) . Class 2 includes objects 
with relatively asymmetric morphology suggesting a post- 
coalescence phase. These systems may have nuclei separated 
by D < 1.5 kpc. 

For those systems with what appears to be a rotating 
disk component, we give in Table [2] their Ha-derived rota- 
tion velocities (uncorrected for the unknown inclinations). 
For all these systems b ut IRAS 03068-5346 and 19297-0406, 
Dasyra et al. ( 2006a b ) measure stellar-derived rotation ve- 
locities from fitting the CO bandheads in their H-band long- 
slit data. Their measurements are systematically lower than 
ours by ^30-50 kms -1 , suggesting that the kinematics of 
the stellar and gaseous components of these rotating disks 
may actually be de-coupled from one another. This de- 
serves further investigation, particularly with spatially re- 
solved stellar velocities at a comparable resolution to our 
gas measurements. 

Assuming virial equilibrium, our gas rotational veloci- 
ties allow us to estimate a mass (M v i r i a i = (2Rv 2 )/G), which 
is given in the following column in Table [2] This, in turn, 
allows us to calculate the Toomre Q ( |Toomre|1964| param- 
eter for characterising the stability of a disk supported by 
differential rotation and random motion against local ax- 
isymmetric perturbations. (A review of the criterion applied 
to gravitationally coupled stars and gas in a disk is given 



in | Jog|p~996[ see also |Downes &; Solomon|[i998 ). This crite- 
rion says whether a disk is stable to collapse (Q > 1) or not 
(Q < 1), and is calculated as following: 



Q- 



nGZ n 



(2) 



where <j r is the uniformly- weighted velocity dispersion of the 
galaxy ( |Davies et al.|201l| , k = 1.5v rot /R, and S maS s is the 
surface mass density (= mass/surface area, where the mass 
is estimated from the Virial mass, M v i r i a i = (2Rv 2 ) / G) . Our 
results are discussed in Section |4~T1 



In Tabled we indicate whether evidence for warm ion- 
ized (wi), neutral atomic (na) or molecular (mol) gas out- 
flows in each galaxy are known either from the literature 
or from this study. A discussion of the differences between 
the outflows measured in the different phases is presented 
in Section 14.31 The next two columns in Table |3] show the 
average centroid velocity (relative to the systemic veloc- 
ity) and line width (FWHM) of the outflow components 
(typical errors of ±30%). The final two columns show the 
two commonly used ways of calculating the "maximum" 
wind speed: Av max ,i = | (v) hroad - Av broad(FWH M)/2| (e.g. 



Veilleux et al. 2005 Rupke et al. 2005b|, and Ai> max ,2 = 


\{v) broad - 2(7broad (e.g. 


Martin 


2005 |Genzel et al.||2011|) 



where we use v 



4 DISCUSSION 

4.1 Gas dynamics and the Toomre Q parameter 

One of the primary goals of studying local ULIRGs is to 
gain insights into the physical processes at work in high-z 
intensely star forming galaxies. What determines the struc- 
ture and nature of the star formation? In the merger pro- 
cess, does the gas become immediately unstable and col- 
lapse in clumps to form stars, or does it get driven to the 
centre to fuel a nuclear starburst and/or AGN? What gov- 
erns these processes and their timescales? For example, con- 
trary to what was expected, the rest-frame UV/optical mor- 
phologies of z > 1 "normal" star-forming galaxies have often 
been found to be dominated by several giant (kpc size) star- 



1996 


Elmegreen et al.|2007||2009; Elmegreen & Elmegreen 


2005 


2006; 


Forster Schreiber et al. 


2009, 


2011 


), implying 



that the gas in these disks forms stars in-situ, rather than 
being driven into the centre to fuel a nuclear starburst. This 
is of interest here as these clumpy, asymmetric structures 
often resemble local (z ~ 0) mergers ( |Conselice|[2003T 



et al. 12004b 



Lotz 



One way of probing what physical mechanisms are at 
work is to measure the stability of a system's disk to gravita- 
tional collapse via the Toomre Q parameter ( Toomre|T9 64 ) . 
In order for star formation to commence, the ISM within a 
rotating disk must become globally gravitationally unstable 
whilst under the influence of shear resulting from differen- 
tial rotation. Regions that become too big will be torn apart 



by shear before gravitational collapse can occur. Elme green] 
( |1999| ) showed that large-scale gravitational collapse must 
set in whenever the gas density becomes a significant frac- 
tion of the "virial" density, so it follows that large, turbulent, 
gas-rich disks are prone to developing star forming clumps. 
SMGs (IChapman et al.||2003|, for example, with SFRs of 



many 100s M yr , have been shown to have Q <C 1 (Tac- 



coni et aL||2006 ). This is unsurprising given that they are 
known to fulfil the criteria for being "maximal starbursts" 
(Tac coni et al.|2006 ). The situation in local ULIRGs is, how- 
ever, less polarised. Q has been found to vary from system 



to system with some <1 and some >1 (Downes & Solomon 



1998) 



In our sample, 3/8 systems with measurable Qs (i.e. 
that exhibit a gaseous disk) have Q < 1 (IRAS 09111- 
1007W , IRAS 20046-0623 and IRAS 20414-1651), whereas 
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Table 2. Galaxy properties of the ULIRG sample. 



Galaxy 


Morph 


Separation of 


^rot 6 


Virial 


Toomre 




classify 


nuclei (kpc) 


(kms- 1 ) 


mass (Mq) 


Q 


IRAS 00198-7926 


1 


6 








IRAS 00335-2732 


2 










IRAS 03068-5346 


2 




60 


6.0 x 10 9 


2.08 


IRAS 05189-2524 













IRAS 06035-7102 


1 


6 








IRAS 09111-1007W 


1 


40 


110 


3.0 x 10 10 


0.93 


IRAS 12112+0305 


1 


4 








IRAS 14348-1447 


1 


4.8 C 








IRAS 14378-3651 


2 










TR AS 1 7208-0014 


2 




140 


3.3 x 10 10 


1.32 


IRAS 19254-7245 


1 


10 








IRAS 19297-0406 


2 




175 


7.1 x 10 10 


1.32 


IRAS 20046-0623 


2 




185 


1.3 x 10 11 


0.87 


IRAS 20414-1651 


2 




150 


5.3 x 10 10 


0.79 


IRAS 20551-4250 


2 




90 


5.5 x 10 10 


1.89 


IRAS 21504-0628 







175 


1.4 x 10 11 


1.01 


IRAS 22491-1808 


2 


1.2, 2.4 d 








IRAS 23128-5919 


1 


4 









a Class 0: single isolated galaxies with no obvious evidence for strong past or ongoing interaction. Class 1: pre-coalescence phase with 

two nuclei. Class 2: relatively asymmetric morphology, or nuclei separated by <1.5 kpc, suggesting a post-coalescence phase. 

b Disk rotation velocity (for objects with clear disks and rotation) 

c |Carico et ah] ( fl990l>; |Scoville| Jl999t; iSurace et al.| (|2000t 

d Triple nuclei flCarico et al.|1990| |Surace et al.|2000| |Cui et al.|200l| ) 

Table 3. Outflow properties of the ULIRG sample. 



Galaxy 


Known evidence 


^outflow k 


^broad(FWHM) C 




Av m ax,2 d 




for out flows a 


(kms" 1 ) 


(kms" 1 ) 


(kms" 1 ) 


(kms" 1 ) 


IRAS 00198-7926 


wi 


-400 


2000 


-1400 


-2100 


IRAS 00335-2732 












IRAS 03068-5346 












IRAS 05189-2524 


wi, mol 


-3000 


4000 


-5000 


-6400 


IRAS 06035-7102 


wi 


-550 


2500 


-1800 


-2700 


IRAS 09111-1007W 


wi 


-250 


700 


-600 


-850 


IRAS 12112+0305 


wi, mol 


-300 


500 


-550 


-725 


IRAS 14348-1447 


wi, mol 


-150 


750 


-525 


-790 


IRAS 14378-3651 


wi, mol 


-300 


700 


-650 


-900 


IRAS 17208-0014 


na, mol 










IRAS 19254-7245 


wi 


>500 


2000 


>-1500 


>-2200 


IRAS 19297-0406 


wi, na, mol 


-0 


800 


-400 


-680 


IRAS 20046-0623 


na 










IRAS 20414-1651 


na 










IRAS 20551-4250 


wi, mol 


-150 


700 


-500 


-750 


IRAS 21504-0628 












IRAS 22491-1808 












IRAS 23128-5919 


wi 


>-450 


700 


>-800 


>-1050 



fl wi = war m ionized (ref: this paper; iLfpari et al . 2003 S poon & Holt||2009||Dasyra et al.||2011| for IRAS 05189-2524); na = neutral 
atomic (ref: [Martin 2005 Rupke et al. 2005b); mol = molecular ( |Sturm et al.|2011| E. Sturm, priv. comm.). 

b Outflow — (^)broad = average centroid velocity (relative to systemic velocity) of the outflow components. Typical errors are ±30%. 
c Average FWHM of the outflow components. Typical errors are ±30%. 

d The "maximum" wind speed calculated in the two ways used in the literature: Av maX) i — |(^)broad — A^bro ad(FWHM)/^| 
|Veilleux et al.|2005| [Rupke et al.|2005b| >, and A^ maX)2 = |(^)broad - 2cr broad | (e.g. |Martin|2005| |Genzel et al.|20lT] >. 
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the remaining five appear to have gravitationally stable disks 
(Q > 1). We note here that our measurements assume that 
the stars and gas are gravitationally coupled (Jog 1996 



Downes &; Solomon|1998 ). This may not be the case (in fact, 



given the disturbed nature of these systems, is highly un- 
likely), however, we do not have measurements of the stellar 
kinematics to quantify what effect this assumption has (the 
Call stellar absorption lines fall just outside of our wave- 
length range at the redshift of our sources). 

This variation in Qs can be interpreted in terms of the 
concept of self-regulated star formation (e.g., Gammie|2001 



Thompson et al.|2005 Tacconi et al.|2006| ). Gravitationally 



unstable disks with Q < 1 collapse and start forming stars 
vigorously. The resulting feedback increases the local veloc- 
ity dispersion which in turn increases Q and has a negative 
effect on the star formation rate. Like this, Q self-regulates 
at Q ~ 1. In each individual case we are likely seeing the 
system simply at a particular point in this self-regulation 
cycle. 

Using detailed IFU observations of star forming galaxies 
at z ~ 1, Genzel et al. (2011) mapped out the Q parame- 



ter on a spaxel-by-spaxel basis and found that giant star 
forming clumps are found where locally Q < 1, even if the 
global average Q > 1. While it is outside the scope of this 
paper to derive spatially resolved Q maps for our targets, 
it is very likely that even systems that we measure to have 
Q > 1 will have locally unstable regions that are forming 
stars. For example, IRAS 20551 has a disk Q parameter of 
2.08, but strong Ha emission is observed in tidal features 
implying strong extra- nuclear star formation. 

The emerging view is that merger-induced starbursts 
have a nuclear component, and also an extended component 
in the form of massive star forming clumps in the disk ( Bour- 
naud|2011| ). What we see in any particular case all depends 



on the exact merger stage and individual characteristics of 
the system. 

4.2 Line ratios in the nuclear and outer regions 

Forbidden-to-recombination line ratios can be useful in di- 
agnosing the strength of the ionizing field strength. In 
this study we have examined the [Nil] /Ha, [Sn]/Ha and 
[Oi]/Ha ratios, chosen to be close in wavelength so as to 
reduce the effects of reddening as much as possible. [S II] 
and [O i] are produced in the partially ionized zone at the 
edge of an ionized nebula. This zone is large and extended 
for hard radiation fields so the [S n]/Ha and [O i]/Ha ratios 
become stronger for Seyferts compared to SB/LINER galax- 
ies (Yuan et al. 2010), or in the presence of shocks ( |Allen 
et al.||2008| ). [Sll]/Ha is particularly sensitive to shock ion- 



ization because relatively cool, high-density regions form be- 
hind shock fronts which emit strongly in [S II] (|Dopita 1997 



|Oey et al.|2000 ). The [Nil] /Ha ratio is more sensitive to the 



presence of a low-level AGN than the other ratios because 
[N n]/Ha is a linear function of nebular metallicity until high 
metallicities where the ratio plateaus. Beyond this, an AGN 
contribution moves the line ratio higher ( Yuan et al.|2010 ). 
In general, the highest line ratios in all three combinations 
are found in the faintest regions, on the outskirts of the 
galaxies, and depressed line ratios in the nuclear regions is 
commonly observed. 

7/18 systems in our sample show a nuclear depression in 



all three line ratios ([Nil] /Ha, [Sll]/Ha and [Ol]/Ha), and 
a further 7/18 exhibit a nuclear depression just in [Sll]/Ha 
and [Oi]/Ha. IRAS 06035, 12112, 19254 and 20046 do not 
show a correlation. What is the cause of this? The existence 
of cases where only two of the line ratios are low excludes 
effects due to Ha stellar absorption artificially reducing the 
ratios, and artefacts due to multi-component fitting, since 
Ha is part of all three ratios. Furthermore, differential ex- 
tinction across the face of the galaxies should also not affect 
the ratios as they are specifically chosen to be reddening free 
( jBaldwin et al.|1981|[Veilleux fe Osterbrock|1987| . Thus we 
must be seeing a real physical effect due to either ionization 
gradients or abundance changes (total metallicity and/or 
relative abundances) between the nuclear and outer parts of 
these systems. 

In order to investigate this phenomenon further, we 
have plotted all three line ratios for CI and C2 in every 
spaxel in each system that shows a nuclear depression (i.e. 
excluding IRAS 06035, 12112, 19254 and 20046) vs. Ha flux 
in Fig. [5] Here we use the Ha flux as a proxy for distance, 
since this avoids any problems due to absolute differences in 
"distance" from system-to-system. Generally Ha peaks at 
the centre and falls of with radius. To aid the eye, we have 
also plotted the sigma-clipped mean ratio and its ±la error 
in bins of 0.15. A clear correlation can be seen between the 
[S 11] /Ha and [O l]/Ha ratios and Ha flux (distance from the 
nucleus), and a weaker one in [Nil] /Ha. The black and red 
points/lines differentiate those systems with log(LiR)<12.1 
and >12.1. Although no significant difference can be seen 
between the two Lir groups in [Sn]/Ha and [Ol]/Ha, this 
shows that systems with depressions in all three ratios (in- 
cluding [Nil] /Ha) tend to have higher Lir. 

All four ULIRGs that do not show any nuclear line 
ratio depressions are interacting double-nucleus or highly 
disturbed systems of class 1 (pre-coalescence phase). Those 
that show depressions in all three ratios are predominantly of 
class 2 (asymmetric morphology and/or nuclei separated by 
D < 1.5 kpc, suggesting a post-coalescence phase). For the 
group of seven that show depressed nuclei in all three ratios, 
five are centrally concentrated systems. There is no partic- 
ular correlation with spectral classification (SB, LINER or 
Sy2). 

In summary, nuclear line ratio depressions appear to be 
found in the less disturbed systems (either isolated class or 
post-coalescence class 2), and a depression in all three ratios 
is associated with high IR luminosity and a central nuclear 
concentration. This implies that the nuclear-to-outskirts line 
ratio changes are predominantly an ionization effect. If the 
system is centrally concentrated then it is likely that the 
nuclear regions are dense and very dusty (therefore highly 
extinguished) , so in the optical we see only the surface layers 
of the gas. Thus one possibility is that the low line ratios 
reflect the fact that this gas is not experiencing the full ra- 
diation field of the starburst/AGN. Inclination effects may 
also play a role, since the nucleus may appear obscured along 
our line-of-sight, but not to gas located further out in the 
plane of the sky. It would be illuminating to repeat this ex- 
periment with near/mid-IR observations where the effects 
of obscuration are mitigated. 

Instead of asking why the nuclear regions have low line 
ratios, one might instead ask why the outer regions exhibit 
higher line ratios. The diagnostic diagrams of Fig. [4a| show 
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that shocks may play a strong part in the ionization of the 
outer (fainter) regions of very many of our sample galaxies. 
Radially increasing line ratios may therefore also result from 
the dilution of the nuclear (SB or AGN) radiation field with 
radius and/or an increasing effect of shock ionization. 

In their spatially-resolved study of LIRG galaxies, 
Monreal-Ibero et al.| (2010) found a correlation between the 
gas excitation (as measured by the [Nil] /Ha, [Sll]/Ha and 
[Ol]/Ha line ratios) and the velocity dispersion of the emis- 



sion lines in interacting/merging systems (see also |Monreal- 
|Ibero et al.|2006 ), and interpreted this relationship as arising 



from an increasing importance of shocks as ionizing sources 
in the extranuclear regions, since shocks can be a cause of 
increased line widths. We, however, find no correlation be- 
tween the three line ratios and the Ha FWHM in our ULIRG 
sample. 

It is not altogether unsurprising that we do not see 
a correlation of this kind (see discussion above). We often 
see the highest line widths in the nuclear regions, precisely 
where the line ratios have a tendency to be at their lowest. 
Given the lack of a clear correlation between the line ratios 
and Ha FWHM, the LIRG relationship found by |Monreal-| 
Ibero et al. (2010) does not straightforwardly seem to apply 



to ULIRGs. 

The other explanation for the line ratio gradients is 
changes in metal abundance, although we have no real means 
to test this with our dataset. Compared to galaxies of the 
same luminosity and mass, LIRGs and ULIRGs are known 



to have an oxygen under-abundance in their nuclei (JRupke 
et al.|2008 ). This is attributed to the merger- induced radial 



inflow of relatively low metallicity gas from large radii to 
the centre, thus diluting the abundance of the nuclear gas 
( Rupke et aT1|2008| [Soto fe Martin||2010| ). A decreased nu- 
clear metallicity would indeed reduce the [N n]/Ha, [S n]/Ha 
and [Ol]/Ha emission line ratios depending on the relative 
abundance of the infalling gas. However, it is difficult to rec- 
oncile the fact that of the three ratios [N n]/Ha is most sensi- 
tive to metallicity yet is the one with the weakest correlation 
with radius, with an abundance explanation. Furthermore, 
if unstable disks (Q<1) are prone to forming massive in-situ 
star forming clumps rather than funnelling most of this gas 
to the centre, then systems with Q<1 might be expected not 
to show depressed nuclear line ratios. This is not the case; 
for example the disk in IRAS 20414-1651 has Q=0.79 (the 
lowest in our sample), but very clearly depressed line ratios 
in both [Sll]/Ha and [Oi]/Ha. 

Thus we conclude that the line ratio gradients are due 
to ionization differences, and very likely reflect obscuration 
in the central regions mitigating the full effects of the star- 
burst /AGN radiation field and/or the increasing dominance 
of shock ionization with radius. 



4.3 Ionized gas outflows 

11/18 of our galaxies show evidence for outflowing warm 
ionized gas with Av max between 500 and a few 1000 kins" 1 . 
Of the SB galaxies in our sample with outflow detections, 
Ai> ma x ranges between 400-800 kms" 1 (depending on the 
method used to calculate Av max ), and some classified as 
Sy2 (i.e. contain an AGN) have Av max > 1000 kms -1 . The 
fact that our Sy2 ULIRGs exhibit faster outflows than the 
pure starburst systems is consistent with previous studies 



that suggest AGN are required to drive winds to ve locities 
>500 kms" 1 jLipari et al.|2003| |Rupke et al.|2005a ). 

Systems that exhibit line components with widths 
>1500 kms -1 include IRAS 00198-7926, 05189-2524, 06035- 
7102, and 19254-7245. In each of these cases the very 
broad components are detected over spatially extended re- 
gions (for more details see the Appendix). Fig. [6] shows the 
Ha+[Nn] line profiles and best-fitting Gaussian models for 
all spaxels in these four systems where one component has 
FWHMM500 kms" 1 . In IRAS 14348-1447 we also find line 
widths >700 kms" 1 . Since all these systems (except 06035- 
7102, but see the Appendix) are known to host an AGN, 
these kind of line widths might suggest we are directly seeing 
evidence of the AGN broad line region (BLR; despite there 
being no type 1 Seyferts in our sample). This cannot be the 
case for two reasons. Firstly, as already mentioned, the broad 
components are detected over spatially extended regions, 
whereas the BLR would be unresolved in our observations 
and therefore its contribution would only be seen in ^one 
spaxel (equivalent to the seeing disk). Secondly, the broad 
components are identified in both the recombination and 
forbidden lines, as can be seen in Fig. [7] where we show si- 
multaneous fits to the [O i] A6300, Ha, [N n] AA6548,6583, and 
[S n]AA6717,6731 lines from two example spaxels in IRAS 
14348-1447 and 19254-7245. The case is the same for the 
other systems. This indicates that any BLR contribution 
must be small since forbidden line emission in the BLR is 
collisionally de-excited. These line widths must therefore re- 
flect very high velocity gas in the narrow line region or halo 
of the galaxy, and is a very important result of this study. 

In Fig. [8] we place our outflow measurements in con- 
text with others from the literature by plotting them as 
a function of Lir. We include data on ULIRGs compiled 
from Nai absorption line (Rupke et al. 2005b a Martin 



2005), mid-IR [Nem] emission line (Spoon fe Holt 2009), 



and far-IR OH 79 |^m ( |Sturm et al. 



2011) measurements, 



and Nai absorption line measurements of LIRGs (Rupke 
et al.|2005b ) and local starburst galaxies (Schwartz & Mar- 



tin 



2004). Where possible, we have recalculated the mea- 



surements given in these papers to correspond to our def- 



inition of Av max ,i (see Section 3.1). This plot shows the 
now well-known trend of increasing outflow velocity with 
IR luminosity ( |Rupke et aLl|2005b| |Martin|[2005l ), despite 
the inhomogeneous tracers used. The outflow speeds of the 
starbursts in our sample are consistent with those found 
by other authors, except for IRAS 06035-7102 whose out- 
flow speed is closer to the Sy2 measurements, suggesting an 
AGN contribution. The outflow speeds of our Sy2 systems 
are consistent with the upper end of those found by [Rupke 



et al. (2005a); IRAS 05189-2524 is notable for having an 



outflow speed significantly faster than the average. 

The warm (and cool) phases of winds are thought to be 
driven by (1) the momentum of expanding hot gas heated by 
SNe and stellar winds (ram pressure; Lehnert & Heckman 



1996 Martin 2005 Veilleux et al. 2005), shocks, or X-ray 



Compton heating if an AGN is present, and/or (2) radia- 
tion pressure from the OB stars in the starburst or AGN 



accretion disk on dust grains mixed with the gas (Martin 
Murray et al. 2005| 2011), and/or (3) magnetocen 



2005 



trifugal forces within an AGN ( |Everett||2005| ). Cox et al 



([2004J hypothesised that the collision of two gas-rich galax- 
ies can also produce outflows driven by the shock heating of 
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gas resulting from the transfer of orbital energy rather than 
supernovae. However, since the shock heating should occur 
in the central few kpc region, its effect is likely indistinguish- 
able from that of the SNe without detailed X-ray observa- 
tions, although some evidence for the role of merger-induced 
shock heating has been found in the nearest ULIRG, Arp220. 



Chandra observations of this source (McDowell et al. 2003) 



revealed two distinct X-ray structures with different spatial 
scales; one suggesting the presence of a superwind, and an- 
other much larger one suggesting shock heating from the 
merger process. 

The theory of ram pressure exerted by a supernova 
heated wind predicts a maximum outflow velocity deter- 
mined by the temperature of the momentum-carrying com- 
ponent of the hot wind. In their samples of starburst dom- 
inated ULIRGs, Martin| (| 2005| found values of At>m aXl 2 be- 
tween 300-750 kms - , and |Rupke et ah] ( |2005b| ) found 
(Av max ,i) = 300-400 kms -1 , both of which are consistent 
with a maximum wind speed of order 500 kms -1 . These 
two studies, however, measured the cool neutral phase via 
the Nai absorption line, thus some caution is required since 
a detailed comparison of column-density-weighted mean ab- 
sorption profiles with rms-density- weighted emission profiles 
may be misleading (even if the absorption and emission lines 
were probing the same gas phase). Radiation pressure the- 
ory does not predict such a maximum wind speed, and it 
has been shown that ULIRGs have continuum luminosities 
that are easily sufficient to accelerate significant columns of 
cool/ warm gas to the observed velocities via radiation pres- 
sure ( |Martin|200"5 ). 

It is very likely that the different mechanisms play a 
larger or smaller role in accelerating different material at 
different radii, times and temperature phases, and that if 
we observe a neutral outflow, it may not necessarily follow 
that we also find an ionized gas outflow. For example, evi- 
dence for a molecular (via OH 79 urn; v max = —370 kms -1 ) 
and neutral (Nai D; tWx = —690 kms -1 ) gas outflow has 
been found in IRAS 17208-0014 ( |Martin||2005| |Sturm et al 
2011), but we find no evidence for it in our optical spec- 



troscopy. Neutral gas outflows have also been identified in 
IRAS 20414-1651 with v max = -230 kms -1 (|Rupke et aL 



2005b), for which we do not see an ionized phase counter- 



part in our data. In e.g. IRAS 19297-0406, however, we 
do find an ionized outflow component that is consistent 
with the neutral component found by Martin (2005 with 
Vmax - -700 kms -1 ). 

As a final note, our spatially resolved data show that 
spatial and spectral properties of ULIRGs are clearly com- 
plex and dependent on the interaction/merger state of the 
system. Furthermore, the extinction in ULIRGs is known 
to be high and very patchy (Garcfa-Marin et al. 2009b). 



One must therefore be careful when making a comparison 
with spatially integrated or luminosity-weighted measure- 
ments (e.g. |Martin||2005| |Rupke et al.|[2005b|a| |Spoon 
Holt||2009| |Sturm et al.||2011| ). However, for the case of the 
outflow properties, in the majority of systems in our study 
we find outflow components coincident with the nucleus, 
meaning that integrated slit spectra of the nucleus would 
generally be sufficient. One caveat to this is that in some 
multiple-nucleus systems with even relatively close separa- 
tions between the nuclei, we find that one nucleus is driving 



an outflow and the other not, meaning it would be important 
to know which nucleus a long-slit observation covered. 



5 CONCLUSIONS /SUMMARY 

We have obtained VLT/VIMOS IFU emission- line spec- 
troscopy of a volume limited sample of 18 southern ULIRGs 
selected with z < 0.09 and S < 10. Our sample covers a 
range of ULIRG types starburst dominated, LINER and 
Sy2 — but not the rare Syl class of which there are no 
\ow-z southern hemisphere examples — , distant pairs, close 
pairs, and fully-coalesced systems). Some systems, despite 
their proximity and ULIRG classification, have received very 
little specific attention outside of the more statistical anal- 
yses of the local ULIRG population in the literature. By 
employing an automated Gaussian line fitting program, we 
decomposed the emission line profiles of Ha, [N n] , [S n] , and 
[O i] into individual components, and have presented here for 
each galaxy spatially resolved maps of the Ha line flux and 
kinematics (FWHM and radial velocity), together with line 
ratio maps of [Nil] /Ha, [Sll]/Ha and [Oij/Ha that can be 
used to diagnose the radiation field strength, and electron 
density maps derived from the [S n] doublet ratio. Together 
with these VIMOS data, we also present corresponding HST 
F814W imaging (or DSS imaging if HST data is not avail- 
able) . 

We discuss the results for each galaxy in some detail 
in the Appendix. Our sample contains galaxies with a wide 
range of morphologies and merger states. Only two are clas- 
sified as single isolated galaxies with no obvious evidence 
for strong past or ongoing interaction, seven are classified as 
being in the pre-coalescence phase with two distinct nuclei, 
and nine are classified as having relatively asymmetric mor- 
phologies and/or nuclei separated by <1.5 kpc, suggesting 
they are in a post-coalescence phase. As is expected, we find 
a large range in gas dynamics and excitations. 

Our overall conclusions are summarised in the following: 

• 1 1/18 of our galaxies show evidence for outflowing warm 
ionized gas with Av max between 500 and a few 1000 kms -1 . 
Of the starburst galaxies in our sample with outflow detec- 
tions, Av max ranges between 400-800 kms -1 (depending on 
the method used to calculate Av max )- Some systems classi- 
fied as Sy2 (i.e. contain an AGN) have Av max > 1000 kms -1 . 
These fast outflows cannot be related to the AGN broad line 
region (BLR) since in all cases they are spatially extended 
and seen in both the recombination and forbidden lines. Four 
of the fastest outflows are found in: 

- IRAS 00198-7926: we find line widths >1000 kms -1 
over a ^3 x 3 spaxel region (~3 x 3 kpc) in a region to the 
west of the northern nucleus. These very broad lines are 
redshifted by 300-450 kms -1 with respect to systemic. 
We also detect broad line components (FWHM=400- 
500 kms -1 ) that are blueshifted by a similar amount 
(300-450 kms -1 ), implying the presence of a bipolar wind 
emerging from behind obscuring material associated with 
the northern nucleus. 

- IRAS 05189-2524: in this known Sy2 system we iden- 
tify four spaxels in Ha C3 with line widths >2000 kms -1 , 
which are all blueshifted by few 1000 kms -1 . Blue wings 
to the [S n] doublet line profile are also seen. These re- 
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suits, together with those at other wavelengths from the 
literature, strongly suggest that this emission originates 
in a fast nuclear wind. 

- IRAS 06035-7102: we find line widths >2000 kms -1 
over nine contiguous spaxels (~2.5 x 2.5 kpc) coinci- 
dent with the western nucleus. These are all strongly 
blueshifted. The classification of this system is uncertain, 
however these results suggest that this western nucleus 
contains a hidden AGN that is driving the outflow. 

- IRAS 19254-7245: we find line widths >2000 kms" 1 
over 34 contiguous spaxels in C2 and C3, coincident with 



the southern nucleus (in good agreement with |Bendo et al. 
2009). This nucleus is known to harbour an AGN that is 
presumably driving this outflow. 

The fact that our Sy2 ULIRGs exhibit faster outflows than 
the pure starburst systems is consistent with previous stud- 
ies that suggest A GN are required to drive winds to velo ci- 
ties >500 kms -1 flLfpari et al.||2003[ [Rupke et al.||2005a| ). 

• Our spatially resolved spectroscopy has allowed us to 
map the outflows kinematically and morphologically, and 
in some cases determine for the first time which nucleus 
is driving the wind and examine the interplay between the 
outflow and the host galaxy. 

• Eight systems have clear rotating gaseous disks, and for 
these we measure rotation velocities, derive virial masses, 
and calculate the Toomre Q parameter ( |Toom re 1964). For 
these galaxies, masses fall in the range 6x10 -1.4x10 . 
Three of the eight have Q<1 and 5 have Q>1. We interpret 
this variation in terms of the concept of self-regulated star 
formation where feedback regulates Q to ~1. 

• We find [Nil] /Ha, [S n]/Ha and [O i]/Ha line ratio gra- 
dients in a significant number of systems in our study. Some 
exhibit gradients in just [Sll]/Ha and [Oi]/Ha, and some 
in all three ratios. Depressed line ratios in the nuclear re- 
gions are found in the less disturbed systems (either isolated 
class or post-coalescence class 2), and a depression in all 
three ratios is associated with high IR luminosity and a cen- 
tral nuclear concentration. We discuss a number of possible 
reasons for this, including ionization effects and abundance 
gradients, and conclude that these gradients are most likely 
a result of an increase in the hardness of the radiation field 
responsible for the observed emission lines with radius, most 
likely due to an increasing contribution of shocks. 

• Neutral and/or molecular gas outflows have been iden- 
tified in some galaxies for which we do not see an ionized 
phase counterpart in our data. The reverse also holds: some 
galaxies where we find an ionized gas wind have not been 
found to have a neutral counterpart. Why this might be is a 
matter of on- going research. It is possible that different wind 
driving mechanisms (e.g. ram pressure vs. radiation pres- 
sure) may be playing larger or smaller roles in accelerating 
different material at different radii, times and temperature 
phases. 

• In many cases, we observe strong Ha emission in tidal 
features/tails implying vigorous star formation is occurring 
there. Since it is known that large stellar structures built 



in tidal debris may end up as tidal dwarf galaxies (Weil 



bacher et al. 2000), it is possible that the collisions that 



form ULIRGs may be a significant source of these types of 
objects. 

Because of the stellar and gas dynamics resulting from 



the major merger (which is the root-cause of the ULIRG 
state), and the interaction between the different gas phases 
and feedback- driving mechanisms from the ensuing power- 
ful starbursts and AGN, ULIRGs exhibit extremely com- 
plex kinematics and gas properties. Even in local infrared- 
luminous mergers such as NGC 6240 or NGC 4038/9, the 
complex kinematics of the ionized gas make it difficult to dis- 
tinguish between wind outflows and tidal motions directly 
produced by the interaction. High spatial resolution, multi- 
wavelength observations of the launching sites of winds are 
required to disentangle the complex kinematics produced in 
ULIRGs (e.g. |Rupke fe Veilleux|2011| ). Furthermore, as this 
study has shown, any system when looked at in detail shows 
peculiarities that need detailed study to understand. This 
makes the case for high resolution follow-up observations 
even clearer. 

The following Appendix contains a discussion of each 
ULIRG system individually, with reference to the existing 
literature. 
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Figure 1. High S/N example spectrum from IRAS 00335-2732 showing the full wavelength coverage of the HR Red grism. Wavelengths 
are not corrected to the rest-frame. The emission lines examined in this paper are labelled. 
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Figure 2. Example single and double Gaussian fits to the Ha and [Nll]A6583 lines from a spaxel in IRAS 20046-0623 (the flux ratio 
of [Nll]A6583 to 6548 is fixed at 3:1 by atomic physics so we do not attempt to accurately model the [Nll]A6548 line profile). The % 2 
for both fits is shown, together with a plot of the residuals (in units of a). For the double-Gaussian fit, the velocity difference between 
CI and C2, and the FWHM of each is shown. These plots demonstrate how two components can be needed to fit a line profile, even 
though their velocity separation is small. In this example hints of a faint blue wing to the lines are visible, although the fitting of a third 
component would not be accepted by our x 2 ratio tests. 
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IRAS 00198-7926 




Figure 3a. Maps showing the Ha flux, FWHM and radial velocity, [Nil] /Ha, [Sll]/Ha and [Oi]/Hq line flux ratios and electron 
density distributions for each identified line component in each of the galaxies in our sample (hatched regions represent non-detections). 
The FWHMs are corrected for instrumental broadening, and the radial velocities are relative to v sys . The top-left panel shows the 
corresponding HST F814W image (for IRAS 00198-7926, IRAS 03068-5346, IRAS 20046-0623 and IRAS 21504-0628 no HST imaging 
is available so we instead show the DSS image); the vertical bar shows a linear scale of 5 kpc. The top-right panel shows the summed 
(C1+C2+C3) Ha flux map. The [S il] and [O i] datacubes were smoothed in the spatial directions by a 3x3 rolling summation to increase 
S/N before the emission lines were fit. This also applies to the electron density maps since densities were derived from the summed 
(all components) [Sn] line ratios. Example Ha+[Nll] line profiles from selected spaxels are shown as insets, together with their multi- 
@n2p$a^A0,aMMSiASt©(^d)^^^SJ'ed lines; wavelengths are rest-frame). Nuclear spectra are labelled "n". The bottom inset in each case 
is the integrated spectrum of the whole system. 
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Figure 3b. 
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Figure 3d. The Ha FWHM and radial velocity maps for C3 are separately scaled as shown to accommodate the very broad and 
blueshifted values. 
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Figure 3g. The Ha radial velocity map for C3 is separately scaled as shown to accommodate the larger values. A blow-up of the region 
enclosed by the dashed rectangle is shown with various components labelled (Nn: northern nucleus; Ng: southern nucleus; Rl: Hn region 
l; |Colina et al.|2000| >. The dashed lines on the inset indicate the extent and direction of the outflow from Nn as found with our data. 
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Figure 3h. 
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Figure 3k. The Ha FWHM maps for C2 and C3, and the radial velocity map for C2 are separately scaled as shown to accommodate 
the larger values. 
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Figure 3m. Near the top-left we show the Murp hy et al.| ( |2001| > proposed merger geometry at approximately the right orientation and 
scale compared to the DSS image above and the Ha maps to the right. For reference, the black outlines on the Ha C2 radial velocity 
map show in which spaxels a third component is also identified. The white cross on the Ha C2 FWHM map indicates the spaxel from 
which the line profile fits in Fig. [2] are from. 
© 2012 RAS, MNRAS 000, [Tpl] 
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Figure 3n. 
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Figure 5. [Nil] /Ha, [Sii]/Ha, and [Oi]/Ha flux ratios vs. Ha 
flux for all spaxels in the subsample of systems that show some 
correlation between these measurements (see text for more de- 
tails). The values in each axis have been normalised to the maxi- 
mum and minimum of the range found in each individual system 
for ease of comparison. The black points/lines belong to systems 
with log(Li R )<12.1 and the red to log(Li R )<12.1. The thick lines 
represent the sigma-clipped mean ratios in bins of width 0.15, and 
the filled areas the zblcr on this mean. 



APPENDIX A: A CASE-BY-CASE 
DESCRIPTION OF THE TARGETS 

ULIRGs are a fairly heterogeneous group of galaxy systems, 
bound together by their high IR luminosities. This classifica- 
tion derives from the dust re-processing of strong UV radia- 
tion produced by a large population of young stars and/or an 
AGN by dust and gas. Although ULIRGs are all thought to 
be colliding/merging systems undergoing a transition from 
gas-rich spirals into low-intermediate mass ellipticals (|Col- 



ina et al.|200T]|Genzel et aI.|2001||Tacconi et al.|2002||Naab 
et al.|2006| ), the individual details of each case vary widely in 
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Figure 8. Outflow velocity, A-u ma x vs. Ljr for systems measured 
in this paper (solid symbols), together with data compiled from 
Nal absorption line ( |Rupke et al.|2005b|a| |Martin|2005| >, mid-IR 
[Ne hi] emission line jSpoon fe Holt|2009| ), and far-IR OH 79 urn 
(Sturm et al. 2 Q11| measurements of ULIRGs, and Nal absorp- 
tion line measurements of LIRGs (Rupke et al. 2005b) and local 
starburst galaxies jSchwartz fc Mar tin 2004). Where possible, we 
have recalculated the measurements given in these papers to cor- 
respond to our definition of Ai; max ,i- Note that the absence - 
or presence - of an AGN in IRAS 06035-7102W needs further 
investigation. 



terms of interaction/merger stage, star formation strength, 
AGN activity, etc., meaning that the reasons for their high 
luminosities also vary. In the following we therefore describe 
the VIMOS emission line maps of Figs |3af[3r| for each galaxy 
in turn, discussing each in detail with reference to the liter- 
ature. 



IRAS 00198-7926 

This system, classified as a Sy2 from optical spectroscopy 
(Vad er et al.||1993| , has received very little previous atten- 
tion in the literature and has never been observed by HST. 
Despite the singular appearance of IRAS 00198-7926 in the 
DSS image, the Ha emission in our VIMOS data is con- 
centrated in two regions implying that this is a colliding 
galaxy system (the projection separation between the Ha 
peaks is 4" or 6 kpc). We detect multiple Ha line com- 
ponents (up to three in the two peaks), the kinematics of 
which are complex and difficult to disentangle. Line widths 
of FWHM>500 kms -1 are associated with the whole north- 
ern sub- system, and in a region just to the west of the 
northern nucleus we find line components with widths up 
to 2000 kms -1 (the colours become saturated on the maps 
due to the upper limit of the scale bar) and redshifts of 
300-450 kms -1 with respect to systemic. Fits to the Ha 
and [N n] lines for these spectra are shown in Fig. [6] Taken 
together, this suggests we are seeing a fast (possibly AGN- 
driven) gas flow here. Immediately surrounding this and to 
the east of the northern nucleus, we also detect reasonably 
broad line components (FWHM=400-500 kms -1 ) that are 
blueshifted by a similar amount (300-450 kms -1 compared 
to systemic) . This pattern implies we may be seeing a bipo- 
lar wind emerging from behind obscuring material associ- 
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Figure 3p. 
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Figure 3r. The black ellipse on the Ha C2 radial velocity map indicates where a third component is also observed. 
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Figure 4a. [N il] A6583/Hai vs. [S ll](A6717+A6731)/Ha flux ratio plots. Each point is colour-coded according to its Ha flux (red=bright, 
blue=faint). As a guide, the photoionization models of |Dopita et al.| ( |2006| > and either the shock models of Alle n et al.| ( |2008| or the dusty 
AGN models of |Groves et al.|| 2004) are overplotted (depending on the spectral type classification of the system). In green (1 Zq) and 
blue (2 Zq) are the |Dopita et al.| photoionization models, with solid lines of constant R (and starburst age ranging between 0.5-3 Myr) 
and dashed lines of constant age (and R ranging between to —6). R is defined as the ratio of the mass of the ionizing cluster to 
the pressure of the interstellar medium. Overplotted in red (1 Zq) and orange (2 Zq) are the shock ionization models of | Allen et al.| 
( |2008| ) (without precursor) for n e = 1 cm -3 , with solid lines of constant magnetic parameter (and shock velocities ranging between 
) and dashed lines of constant shock velocity (and magnetic parameter B/n ' rangin 

but are offset by approximately —0.2 dex in both line ratios. On selected 



100-500 kms- 

for shocks with precursors cover a similar area in this diagram, 
plots, instead of the shock models we plot the dusty AGN models from|Groves et aL](|2004| for a range of metallicities (1, 2, and 4 M©, 
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Figure 4b. 
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Figure 4c. 



ated with the northern nucleus. We therefore propose that 
the northern system hosts the Sy2 nucleus, and is deserving 
of more detailed follow-up observations. 

No obvious evidence for disk rotation is apparent in ei- 
ther sub-systems, and there is no radial velocity difference 
between the two nuclei. An arm of Ha emission appears to 
extend from the internuclear region to the east then to the 
north around the northern nucleus. A continuous radial ve- 
locity gradient along this arm implies that it is a coherent 
structure, although it begins blueshifted ^150 kins" 1 com- 
pared to systemic and ends redshifted by ^250 kms -1 . 

The same region surrounding the northern nucleus ex- 
hibiting broad lines has generally high [Nil] /Ha, [Sn]/Ha 
and [O i]/Ha ratios that are consistent with those produced 
by AGN excitation (Fig. 4a) The two nuclei show lower 



[Sn]/Ha and [Ol]/Ha line ratios than the surrounding gas. 
Electron densities peak at ~700 cm -3 in the internuclear re- 
gion, but remain high (>300 cm -3 ) around both the nuclei. 



IRAS 00335-2732 

This starburst dominated galaxy has also received very lit- 
tle previous attention in the literature. In the HST imag- 
ing it appears to be an asymmetric, post-coalescence phase 
system with one prominent tail with a number of conden- 
sations. This tail is not detected in our emission line data. 
However, we do detect emission from a much fainter tail 
seen in the F814W image, extending ~10 kpc to the north- 
east. This Ha-bright tail is redshifted compared to the main 
disk, exhibits low Ha line widths and low line ratios. No 
signature of rotation is seen in the main disk. The [S n]/Ha 
and [Oi]/Ha line ratio maps are depressed in the central 
region of the disk and higher on the outskirts, whereas the 
[Nil] /Ha line ratios are fairly constant across the disk. Gas 
densities of >500 cm -3 are observed in the nuclear region. 



IRAS 03068-5346 

This is another starburst dominated galaxy that has re- 
ceived very little previous attention, and has never been ob- 
served by HST. A clear rotation gradient can be seen within 
the main body of the galaxy implying circular motions of 
vsini — 50 kms -1 (the inclination is unknown). In the DSS 
image, a disk and prominent tail extending 10-15 kpc to the 



south-west can be seen. In our Ha emission-line image this 
south-western tail is also seen, but at a distance of ~10 kpc, 
a bend in the tail towards the east is observed. In the velocity 
map, the tail is a coherent structure that appears attached 
to the galaxy disk and becomes increasingly redshifted to- 
wards the south-west. Such stong Ha emission in the tail 
implies significant SFRs (>few yr _1 ). These high rates 
could be associated with the crossing orbits in very young 
tails (cf . |Toomre fe Toomre|1972l |Wallin||1990| ) . 

A ^5x3 kpc region to the north-west of our IFU 
data exhibits very large Ha line widths (FWHM ~ 500- 
1600 kms -1 ), and is roughly coincident with an area of red- 
shifted gas. All three line ratios are relatively high in this 
whole region. Since this region is on the outskirts of the 
system (>5 kpc from the nucleus), we speculate that here 
we are seeing the impact site of infalling matter. Of course, 
higher resolution observations would be needed to investi- 
gate this further. Multiple line components can be identi- 
fied in the main disk and inner parts of the tail, but none 
are suggestive of outflowing gas. The line diagnostic plot in 
Fig. [4a| suggests that most of the ionized gas in this system 
is excited by young stars, with only a small contribution 
from shocks in the fainter, outer parts, consistent with the 
proposed infall region. 



IRAS 05189-2524 

IRAS 05189-2524 has been classified as a late stage merger 
with a single, compact and very red nucleus and a Sy2 op- 
tical spectrum (Veilleux et al. 1995|. The X-ray spectrum 



dRisaliti et al.|2000| |Imanishi et al.|2003) , SE P fitting jFar"- 
)03|), 



rah et al. |2003) , and Spitzer IR spectroscopy (iVeilleux et al 
2009) all suggest this object contains both a starburst and 
AGN where the AGN provides a significant fraction (^70%; 
Veilleux et al.|2009 ) of the total IR luminosity. The detection 



of N IV] A1488 in the UV by Farrah et al. (2005), however 



implies a large population of W-R stars and that a young 
(3-4 Myr) starburst is certainly present. The HST F814W 
image shows that the compact nucleus of IRAS 05189-2524 
is surrounded by a large system of faint concentric annular 
shells, as is found in ~10% of elliptical galaxies (Schweizer 
1980||Malin fc Carter |1983| . These shells are believed to be 
the remnants of a small galaxy after a head-on collision with 
the larger system |Quinn|1984| |Merrifleld Kuijken|1998 ). 
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Figure 7. Simultaneous fits to the [Oi]A6300, Ha, [N n] AA6548,6583, and [S n] AA6717,6731 lines in the nuclear spectra of IRAS 14348- 
1447 and IRAS 19254-7245, showing how the broad component is present in both the recombination and forbidden lines. (For IRAS 
19254-7245 we can also fit the faint [Ol]A6360 line.) The insets show an enlargement of the [Sll] region. In purple we show the individual 
[Sll] 6717 and 6731 elements of the broad component fit since our routine fits the two lines simultaneously and hence displays it as one 
blended profile. 



No evidence of rotation in seen in either the Ha CI 
or C2 components, implying that we are either seeing a 
disk fully face-on, or that this system is completely pres- 
sure supported. The latter is more likely considering its 
early- type morphology. A relatively broad (500-800 kms -1 ) 
and blueshifted (^350-400 kms -1 with respect to CI) 
Ha C2 line component is identified over a ^5 kpc 2 circu- 
lar region centred on the nucleus. These velocities are in 
good agreement with the blueshifted high-ionization lines of 
[O ml, [Ariii]AA7137,775 2, and [S m]AA9069,9531 measured 
by Rupke et al. (2005a) to have velocities of —500 kms -1 
with respect to systemic, and the similarly blueshifted 
(—400 kms -1 ) Nai absorption. The very large area over 
which we see this broad, blueshifted emission indicates a 
significant wind outflow with a wide opening angle that is 
perhaps directed along our line of sight. 



therefore be coming from t he wind) 2 ] Tog ether, the results 
oflSpoon & Holtl (12009} andlDasyra et al.l 12011) imply that 



Rupke et al. (2005a) noted that some of the high ion- 



ization lines also have even broader blue wings. We also find 
five spaxels in the very centre exhibiting very broad wings 
to the Ha, [Nil] and [Sll] emission lines. Here we fit a third 
component with widths >1500 kms -1 and centroid radial 
velocities that are blueshifted by a few 1000 kms -1 (with 
respect to CI). Fits to the Ha and [Nil] lines for these spec- 
tra are shown in Fig. [6] Farrah et al. ( 2005 ) also found broad 
components in R/3 and [O III] A5007 in their HST/STIS spec- 
troscopy of the nucleus, with systematic blueshifts compared 
to the other optical lines of 300-1200 kms -1 . In their Spitzer 
spectroscopy of this system, Spoon & Holt ( 2009| ) identified 
strongly blueshifted [Nev] emission (^1000-1500 kms -1 ). 
Dasyra et al. (2011) also detected blueshifted [O iv] with 
even higher velocities. |Rupke et al.| ( [2005a ) find a very good 
match between their Nai absorption profile and the [Nil] 
and [Oiii] emission line profiles (all showing blue wings), 
thus providing very strong support that these features origi- 
nate in a wind. Strongly blueshifted molecular line emission 
is also seen in Hers chel /P ACS spectroscopy of this system 
(E. Sturm, priv. comm.). That these blueshifted wings are 
found in the recombination and forbidden emission lines, 
and absorption lines, rules out an AGN BLR origin. Most of 
the wing emission in the hydrogen recombination lines must 



this wind is accelerating, since higher ionization lines are ex- 
cited nearer to the ionizing source (the ionization potential 
of [Oiv] is 55 eV compared to 97 eV for [Nev]). 

These findings are reminiscent of those of Ru pke <£] 



|Veilleux| (|2011| for Mrk 231, in which they find both a fast 
blueshifted neutral (Nai) outflow extending 2-3 kpc from 
the nucleus in all directions that they associate with the 
AGN, and a slower, starburst wind arising in the extended 
galaxy disk. Thus, IRAS 05189-2524 is an ideal candidate 
for follow-up, spatially resolved observations of the neutral 
and molecular gas phases. 

Finally, the [Nil] /Ha, [Sn]/Ha and [Oi]/Ha line ratios 
are all very high in the outer disk (see Fig. |4a[ ). These are 
in fact the highest line ratios found in our whole sample, 
and are outside the range predicted by the photoionization, 
shock or AGN models. These results may reflect either a 
particularly strong influence of the AGN ionization on the 
whole system, or an unusually high metallicity (>4Z©) All 
three line ratios are much lower in the very nuclear region, 
but still remain above those predicted by just photoioniza- 
tion. 



IRAS 06035-7102 

IRAS 06035-7102 is a double-nucleus system. The separation 
of the nuclei (6 kpc) and the strong tidal features suggest 
that this system is very close to coalescence. This is sup- 
ported by the smooth Ha radial velocity map, indicating 
the two colliding disks are in the process of forming a single 
structure. From SED fitting, Farrah et al. (2003) classified 



it as a SB/Hn system, however |Weedman &; Houck (2009) 
included this source in their compilation of obscured AGNs 
due to the strength of its silicate absorption and weakness of 



2 Results from previous studies ( |Young et al.|l996||Veilleux et ah] 
|1999| |Farrah et "aT]|2005) imply that the BLR is moderately, 
although not completely, obscured, meaning that some of the 
broad recombination line (Ha, H/3) emission in the central nu- 
clear spaxel may originate here. 
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its PAH emission found from Spitzer /IRS spectroscopy. Our 
spatially resolved kinematics show that the western nucleus 
exhibits very broad, blueshifted wings to the Ha line pro- 
file that are suggestive of the type of fast outflows usually 
associated with AGN. The third Ha component, identified 
in nine contiguous spaxels (^2.5x2.5 kpc) surrounding the 
western nucleus, has line widths >1000 kms -1 and centroid 
velocities that are blueshifted by 500-600 kms -1 with re- 
spect to CI. Those seven with line widths >1500 kms -1 
are shown in Fig. [6] Using the outflow speed formulation 
described in Section |3.1| these measurements convert to 
= 1800 kms . From Fig. ^ we see that this 
value puts IRAS 06035-7102 significantly outside the range 
of outflow velocities measured for other starburst dominated 
ULIRGs. We therefore conclude that the most likely situa- 
tion is that the outflow from IRAS 06035W is being driven 
by a heavily obscured AGN, although we cannot be sure 
without further observations to investigate the presence of 
an AGN. 

Broad Ha CI line widths (200-450 kms -1 ) are found 
in the fainter regions just to the north of both nuclei be- 
tween the tidal tails, and are coincident with elevated line 
ratios in [Nil] /Ha, [Sn]/Ha (only one component of [S n] 
could be identified due to strong sky line residuals at this 
wavelength). Most of these Ha- faint regions have line ratios 
that fall in the shocked region of diagnostic diagram shown 
in Fig. |4a| Unlike the majority of the system where low gas 
densities are measured, large densities (few 1000 cm -3 ) are 
found in the north-eastern tidal arm. 



IRAS 09111-1007W (SB+Sy2) 

IRAS 09111-1007 consists of two widely spaced but interact- 
ing galaxies. The two nuclei have a projected se paration of 
~40 kpc and a velocity difference of 425 kms -1 (Due et al. 



1997). Monte Carlo simulations give the probability of the 



pair being bound as 0.88 ( jSchweizer 1987 ). The western com- 
ponent dominates the far-IR flux at both 60 |^m and 350 |^m, 
carrying ^80% of the total system luminosity (|Khan et al 



2005). The wide separation and the IR luminosity-to-H2 
mass ratio suggest that the pair are possibly at an early 



stage of interaction (Khan et al. 2005), although the high 



IR luminosity of the system might suggest otherwise. Com- 
pared to the eastern galaxy, the western one seems to have 
a high dust content since H/3 is invisible and there is strong 
Nai absorption ( |Duc et al~|l997 ). 

Observations of both galaxies in the system were ob- 
tained, but here we only discuss the western galaxy since 
the IRAS 09111-1007E is not classified as a ULIRG. IRAS 
09111-1007W is a ULIRG with a SB+Sy2 spectrum. A study 
of the two together will form the subject of a future paper 
(Khan et al., in prep). 

The HST F814W image shows that IRAS 09111-1007W 
is a barred spiral. Our Ha radial velocity map shows a clear 
disk rotation pattern, and distortions in the central region 
follow the bar shape of the bar. A second velocity compo- 
nent is observed in the central ^5 kpc, and it is likely that 
this component is kinematically related to the bar. Ha line 
widths across the whole disk are low (<200 kms -1 ), and 
the distribution of line ratios plotted in Fig. [4a] imply that 
this galaxy has a high metallicity. In CI, the [Nil] /Ha line 
ratio remains low across the whole disk. This is in contrast 



to the [Sii]/Ha ratio, which shows a strong radial gradi- 
ent increasing at larger galactocentric distances (which is to 
some extent mirrored in the [Oi]/Ha ratio). A few spax- 
els covering the nucleus exhibit broad (^500 kms -1 ) and 
blueshifted (200-300 kms -1 with respect to CI) emission, 
which we have assigned to C3. Since these components also 
have very high [Nil] /Ha ratios (but are undetected in [S II] 
or [Oi]), we conclude that this component represents emis- 
sion associated with a nuclear wind outflow. The velocities 
would imply a starburst-driven wind rather than an AGN 
outflow. 

IRAS 12112+0305 



IRAS 12112+0305 is classified as a LINER (Veilleux et al 



1999) , and through SED modelling [Far rah et al.|fl2003| ) were 
able to say confidently that this object is powered by a 
starburst with no evidence for an AGN component. Mid-IR 
spectroscopy from Spitzer suggests that the AGN fraction 
is 00% jVeilleux et al.|2009) . IRAS 12112+0305 was pre- 
viously studied with integral field spectroscopy by |Colina| 



etaL] ( pOOOf using INTEGRAL on the WHT. They found 
that the observed ionized gas distribution is decoupled from 
the stellar main body of the galaxy. They also confirmed that 
the optical properties of the two dust enshrouded {Ay — 3- 
8 mag) nuclei are consistent with being weak [O i] LINERs, 
and measured the integrated star formation rate of the star- 
bursts associated with the two nuclei to be ^80 Mq yr -1 
(this is significantly less than the IR-derived value given in 
Table [I] perhaps reflecting that the dust is optically thick 
at Ha; e.g. | Rang wala et al.||2011| ). 

Our significantly higher spatial resolution data clearly 
resolve the two nuclear complexes, the main Ha-bright H II 
regions, and the southern arc-like tidal tail. The Ha velocity 
field is quite distorted across the main body of the galaxy. 
The gas in the southern tidal tail is quite well ordered, and 
becomes progressively and smoothly blueshifted towards the 
south-east, showing that this tail is indeed attached to the 
central regions. It's Ha brightness would suggest that stars 
are forming vigorously here. The patchy observed density 
distribution shows that some areas outside the main nuclear 
collision region and Hn complex have quite high densities. 
The young massive Hn region Rl (Colina et al. 2000[ la- 



belled on the HST image) appears to be forming within a 
quiescent region of the collision since it is located in a region 
of narrow line widths, low [Nil] /Ha, [Sn]/Ha and [Oi]/Ha 
line ratios, and a comparatively low gas density. In contrast, 
the southern nucleus is embedded in a region of broader line 
widths (300-400 kms -1 ) and higher line ratios. 

Redshifted (200-400 kms -1 relative to v sys ) and broad 
(500-600 kms -1 ) Ha emission is seen in spaxels extending 
to the north-east from the northern nucleus (Nn; shown on 
the HST image with dashed lines). This suggests we are 
seeing an outflow from the northern nucleus and H n region 
complex. Herschel /PACS spectroscopy also reveals evidence 
for outflowing molecular gas in the OH lines (E. Sturm, priv. 
comm.). 



IRAS 14348-1447 

IRAS 14348-1447 consists of two nuclei separated by 4.8 kpc 
dCarico et aT][l990) |Scoville|p99| |Surace et aH|205o) with 
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a bright tidal tail to the north and a fainter one to the 
south-west ( |Sanders et"aT]|1988| [Melnick fe Mirabel||1990 ). 
According to previous work, this system is one of the most 
luminous and molecular gas rich ULIRGs that show no clear 
evidence of AGN activity (i.e. broad line emission or strong 
high ionization lines; Veilleux et al.|1 997; Eva ns et al.|2 00Q). 
However, SE P fitting (|Farrah et al.||2Q03| ) and hard X-ray 
observations ( Franceschini et al.||2003] fTmpry that this ob- 
ject contains an AGN, but the total luminosity is dominated 
by the starburst (the AGN contribution is estimated to be 
-17% of the total IR flux; |Farrah et al.|2003| |Veilleux et all 
2009). Indeed, the south-western nucleus has been classified 



as a Sy2 b y some authors (Sanders et al. 1988 Nakajima 
et al.fl991 ). Our Ha maps show that broad (~700 kms - ) 
and blueshifted (100-250 kms -1 ) emission in C2 is asso- 
ciated with the south-western nucleus, implying that this 
object is driving a wind outflow. An example spectrum and 
mult i- Gaussian fit from this region is shown in the figure in- 
set, and in Fig. [7] we show a simultaneous multi-component 
fit to all the observed emission lines. The very broad width of 
this blueshifted component suggests that the obscured AGN 
may be the powering source. Indeed, similar molecular-phase 
outflows can be identified in Herschel /P ACS spectroscopy 
(E. Sturm, priv. comm). 

The Ha intensity map shows extended clumpy emission 
throughout the body of the system and tidal tails, in good 
agreement with what was observed by (iMihos & Bothun 



1998). Strong Ha emission is associated both with a num- 



ber of knots to the south-east of the two nuclei at radii 
of (7-20 kpc), and with the northern tail, imply- 

ing that there is vigorous star formation occurring through- 
out the merging system. An interesting arc-shape string of 
Ha clumps to the south seems kinematically distinct (red- 
shifted) from rest of ionized gas in system, in agreement 
with what was found by |Mihos &; Bothun] ( 1998] ). These au- 
thors suggested that these Ha bright clumps may represent 
the early stage of the formation of tidal dwarf galaxies. El- 
evated line ratios and Ha line widths are found in a strip 
between the two nuclei, that are reminiscent of the shock 



regions in Stephan's Quintet (Pietsch et al. 111997 Guillard 



et al. 2009). 



IRAS 14378-3651 

This system appears on the HST F814W imaging morpho- 
logically similar to IRAS 05189-2524, in that it has a single, 
compact nucleus surrounded by a system of faint concen- 
tric shells. It also displays no signature of rotation in the 
Ha kinematics, implying that it is an early-type system. 



Through SED fitting, Farrah et al. ( 2003 ) found no evidence 



for an AGN in this galaxy, although mid-IR line diagnostics 
imply that there is some AGN component ( |Farrah et al. 
2007| |Veilleux et"aI1|2009| |Sturm et aLl|2011| ). A molecular 



gas outflow of 800-1000 kms was observed via the OH 



79 urn and OH 119 U-m emission lines with Herschel (Sturm 
et al.|20lT ), and these authors argued that such high veloci- 
ties implies that an AGN is powering the outflow. In our Ha 
maps we find broad (500-800 kms -1 ) and blueshifted (200- 
400 kms -1 with respect to CI) emission coincident with 
nucleus and elongated in an east- west direction. This com- 
ponent exhibits high [Nil] /Ha ratios but is not elevated in 
in [S n]/Ha or [O i]/Ha. In the line ratio diagnostic diagram 



of Fi g. |4b| we see t hat the highest metallicity AGN models 
from Groves et al. (2004) still do not cover this part of the 
plot, suggesting either metallicities >4 Z©, or some peculiar 
ionization effects. Of the very few systems that have a mea- 
surement of the outflow speeds in both the warm ionized and 
cold molecular phases, this galaxy presents a very unusual 
case where the molecular material is outflowing faster than 
the ionized material. 

In CI, lower line ratios are consistently found in the 
nucleus, whereas much higher values are found on the out- 
skirts. From the diagnostic plot in Fig. |4b[ these values are 
consistent with high metallicities and AGN/shock excita- 
tion. 



IRAS 17208-0014 

This object consists of a single nucleus surrounded by a 
disturbed disk containing several compact Hn regions. Its 
pure SB classification is supported by optical spectroscopy 
dVeilleux et al.|1999|), near-IR imag ing flScoville et al.|2000| ), 
and SED fitting ( Farrah et al.|2003 ), which all argue against 
the presence of an AGN. However, a significant hard X-ray 



source has been detected by both XMM-Newton (Frances- 
Ichini et al.||2003|) and Chandra flPtak et al.||2003|). Mid-IR 



AGN diagnostics predict that an AGN may contribute ^10% 
to the total luminosity ( |Sturm et al.||20i"T ). Our Ha kine- 
matics show an ordered radial velocity field with a clear sig- 
nature of rotation, but distortions that are consistent with 
the disturbed morphology. Ha CI is fairly broad throughout 
the whole disk (200-350 kms -1 ) implying significant turbu- 
lence presumably driven by the intense star formation. A 
second Ha component, seen only in the central regions, is 
much narrower, but not significantly offset in velocity-space. 
IRAS 17208-001 4 was also observed with IFS by |Arribas| 
&; Colina (2003), who also did not detect any evidence for 
an ionized gas outflow. Martin (2005), however, identified 



a neutral gas outfl ow via Nal a bsorption lines with speeds 
>500 kms -1 , and Sturm et al. (2011) measured a molecu- 



lar outflow from Herschel /P ACS observations with speeds 
of few 100 kms -1 . A large area of depressed line ratios con- 
sistently seen in all three diagnostics is coincident with the 
nuclear region. Higher ratios seen in the outer parts of the 



galaxy are consistent with shocks (Fig. 4b). 

Our higher resolution kinematic results support the con- 
clusions of Arribas &; Colina (2003| who suggest that IRAS 



17208-0014 is a luminous, cool ULIRG that is in the final 
coalescence phase. Furthermore, the fact that there is lit- 
tle evidence that an AGN has developed is consistent with 
the scenario put forward in |Colina et al. (2001) that cool 
ULIRGs are the merger of two low-mass (<m*) disk galax- 
ies that will never evolve into a QSO. 



IRAS 19254-7245 (the "Superantennae") 



IRAS 19254-7245 was the subject of study in [Paper I| Here 
we present a brief re-analysis of these data for both the pur- 
poses of completeness, and to concentrate on the evidence 
for outflows. In this study we have re- fit the emission lines 
with the method used for all the other targets to ensure con- 
sistency. The results are, or course, in very good agreement 
with what we found before. The main difference is that with 
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the line- fitting method used here, we were able to identify 
an additional Ha emission line component in some regions 
that is clearly related to a wind outflow. 

The "superantennae" nomenclature refers to two 
extremely long tidal tails, reminiscent of the tails 
in NGC4038/9 (the "antennae galaxies"), that extend 
>500 kpc on either side of the colliding nuclei. The nuclei 
themselves are separated by 8'/ 5 (10 kpc). To form these fea- 
tures a coplanar collision between two massive disk galaxies 
is implied (Melnick &; Mirabel 1990). The southern nucleus 
has been optically classified as a Sy2 and the northern as a 
pure SB (|Mirabel et al.|1991| ). [Berta et al.| ( |2Q03| ) ( [Frances- 1 
chini et al.|2003 ) found that the near-IR (X-ray) luminosity 
from the southern nucleus is dominated by the emission from 
an obscured AGN, and showed that this nucleus is respon- 
sible for about 80% of the total infrared luminosity of the 
system. Both studies also found that the northern nucleus 
does not show evidence for AGN emission and appears to 
be in a post-starburst phase. High angular resolution HST 
imaging shows that a double nucleus may be present in both 
the north and southern components of the Superantennae 
( Borne et al.||1999| , suggesting a multiple merger origin of 
the system. 

In our VIMOS data, the strongest Ha emission is con- 
centrated around the two nuclei, and in particular the south- 
ern Sy2 nucleus. In |Paper l| we clearly identified the rota- 
tion signature of the two disks. The major axis of rotation 
of the southern disk is aligned in a NE-SW direction, and 
the northern disk in a E-W orientation; both disks rotate in 
the same sense (with their northern sides receding; see also 



Mihos & Bothun 1998 ). This is consistent with the aforemen- 



tioned coplanar collision scenario. In |Paper l| we identified 
broad spectral line emission from the core of the southern 
nucleus, and potential outflows with shock-excited spectral 
features near both nuclei. We estimated that <10% of the 
24 urn flux density originates from star formation, imply- 
ing that most of the 24 urn emission originates from the 
AGN in the southern nucleus. In this paper we have re- fit 
the emission lines with the method used for all the other 
targets to ensure consistency. The results are, or course, in 
very good agreement with what we found before. With our 
updated fitting method we identify a third Ha line compo- 
nent in a cluster of spaxels covering the nuclear regions of 
the southern disk. In 8 of these spaxels, two broad compo- 
nents with FWHM of 1500-2500 kms -1 and velocity offsets 
of 600-800 kms -1 are required to fit the asymmetric wings 
of the profile. All those spectra in IRAS 19254-7245 exhibit- 
ing line widths >1500 kms -1 are shown in Fig. [6] and in 
Fig. [7] we show a simultaneous multi-component fit to all 
the observed emission lines in the spectrum of the southern 
nucleus. These line widths are consistent with those found 
in |PaperTl and as suggested there, the presence of these ex- 
tremely broad components and their velocity offsets imply 
the presence of high- velocity outflowing ionized material as- 
sociated with the southern (Sy2) nucleus. Although we mea- 
sure a Av max of 1500-2200 kms -1 (Table^, emission in the 
wings of the Ha- [N n] complex extends from —3300 kms -1 
to +3700 kms -1 at zero intensity. Both [Colina et a~ (1991) 
and Mihos &; Bothun| ( 1998 ) also reported evidence of out- 



flowing material associated with the southern nucleus. 



broadest widths seen across the southern disk and close to 
the nucleus of the northern disk. This implies a very tur- 
bulent ISM caused by the intense feedback. Line ratios are 
generally high, with evidence for ionization from hard radi- 



ation sources (AGN, shocks; Fig. 4b Paper I ) 



IRAS 19297-0406 

HST imaging shows that this object has a complex, com- 
pact nucleus, composed of several luminous cores embedded 
within a common, diffuse envelope (Bushouse et al. 2002). 
One obvious large tidal tail extends to the east away from 
the nucleus. Due to a pointing error, the nucleus of IRAS 
19297-0406 fell on the very edge of our VIMOS field-of-view, 
meaning we only cover the western half of the galaxy (as 
can be seen in Fig. |3l|. Two Ha line components are iden- 
tified in the nuclear regions. Here line widths of the fainter 
component (C2) extend up to FWHM=800-900 kms -1 , but 
there is little evidence for any radial velocity offset between 
the line centroids. Nevertheless, these line widths suggest 
the presence of an outflow, and are surprisingly broad for 
a system without a previously identified AGN component. 
|Martin| ([2005 ) measured a Na I neutral gas outflow from this 
system with speeds >500 kms -1 and Herschel /P 'ACS spec- 
troscopy of the OH 79 |^m line reveals evidence for a slightly 
slower molecular gas outflow (E. Sturm, priv. comm). 

A radial velocity gradient across the system in the 
north-south direction implies that it is in the late stages 
of its merger and is in the process of forming a coher- 
ent disk. The measured rotation velocity is ~175 kms -1 . 
The [Nil] /Ha, [Sn]/Ha and [Oi]/Ha ratios are all signifi- 
cantly depressed in the centre of the disk, with the ratios 
increasing with radius into the shock region of the diagnos- 
tic diagram shown in Fig. |4b| Electron densities throughout 
the main disk are high with values >300 cm -3 , peaking at 
-1000 cm" 3 . 



IRAS 20046-0623 

This optically classified SB system was studied in some de- 



tail my Murphy et al. (2001 ) via optical continuum and Ha 



imaging and Paa integral field spectroscopy. Through anal- 
ysis of its morphology and velocity maps, they proposed a 
merger geometry model, which we have reproduced to scale 
in Fig. |3m| They propose that this is a system of two merging 
disk galaxies: the eastern galaxy is in the background, and it 
is this galaxy to which the northern tidal tail belongs, while 
the foreground western galaxy is observed nearly edge-on 
and could be providing a screen of extinction that obscures 
emission from the central regions of the eastern galaxy. In- 
deed, two nuclei are visible in H-band images of this source 
(Das yra et al.|2006a| . The tidal tail extending to the north 
shows appreciable line emission. From their Ha imaging, 



Murphy et al. (2001 ) estimate star formation rates of a few 



solar masses per year. The relatively short length and high 
surface brightness of the tail support the idea that this sys- 



tem is being viewed soon after its first encounter i Murphy 
et al.pOOip . 

Our VIMOS data are in agreement with the Paa 



Over the rest of the system, (as we found in Paper I) 
the line width of CI varies from —100-700 kms -1 , with the 



data of Murphy et al. (2001), although in our observa- 



tions the Ha redshifted radial velocities in the eastern half 
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of the system are not centrally concentrated, but are dis- 
tributed into two peaks on either side of the disk. One is 
at the base of the northern tidal tail. The Ha line widths 
are FWHM<150 kins" 1 in the tail, but increase above 
200 kms" 1 in the disk and peak at >800 kms" 1 in a re- 
gion to the south of the main disk. This is approximately 
coincident with a peak in all three line ratios. This southern 
region exhibits a strong radial velocity gradient in the east- 
west direction in the same sense as the main disk, implying 
that it is part of the same structure. It is difficult to recon- 
cile the high line widths and their location with the Murphy 
et alT] merger model. 

In the disk and inner region of the northern tail, two 
Ha line components are observed, and in the central part of 
the disk, a third component can be identified. An example 
spectrum and mult i- Gaussian fit from this three-component 
region is shown in the figure inset. An interpretation of 
what is happening in this region is not straightforward and 
would require higher resolution observations to disentangle. 
Components 2 and 3 trace a region of redshifted (C2; 300- 
350 kms" 1 ) and blueshifted (C3; -50 kms" 1 ) emission with 
the velocity gradient oriented in the north-south direction. 
The black outlines on the Ha C2 radial velocity map indi- 
cate where the third component is identified. At the same 
location, the CI component is at its broadest. Taken to- 
gether, this may suggest the presence of an outflow with 
velocities of ^200 kms" 1 , which is is consistent with the 
Nai outflow found by Rup ke et al.| ( 2005b| with a velocity 
of ^220 kms" 1 . More detailed, higher resolution data are 
needed to disentangle the complex kinematics of the merger 
and confirm the presence of an outflow. 



IRAS 20414-1651 

This system has been imaged on several occasions from the 



ground, from optical to near -IR wavelengths ([Melnick 
Mirabel|1990||Duc et al.|1997| [Surace et al.|2000| |Bushouse 



et al.||2002[ ). HST images reveal that this system is com- 
posed of a single central body that is highly elongated, with 
hints of faint, amorphous emission surrounding it in various 
places. From these previous data it is unclear whether the 
south-western component, separated from the main galaxy 
body by 10'/ 5 (~18 kpc) is tidally connected or related to 
the main galaxy disk. Our observations clearly detect this 
companion in Ha (implying strong star formation is taking 
place) , and since the radial velocity of the emission line gas is 
redshifted by only ^300 kms" 1 compared to the main disk, 
this suggests that the two are connected and interacting. In 
the I-band the central body has a complex morphology, and 
no clear sign of any dominant nucleus. In the H-band, how- 
ever, the morphology is quite regular, showing a distinct cen- 
tral core (Bus house et al.|2002 ). These near-IR observations 
imply that dust extinction in the main disk is high, and are 
consistent with the optical starburst classification. Indeed, 
four out of five mid-IR indicators as seen with Spitzer are 



consistent with there being no AGN contribution (Veilleux 
et al.||2009| ) 



In Ha a clear signature of rotation is seen across main 
disk. The [S n]/Ha and [O l]/Ha ratios are both significantly 
depressed in the centre of the disk, but rise into the shock 



regime (Fig. 4b) in the outer parts. This is not mirrored 



in the [Nil] /Ha ratio, which remains fairly high across the 



whole disk. The line ratios in the south-western companion 
are all low, whereas the gas densities are much higher than 
in the main galaxy disk. 

A region of broad Ha emission (250-450 kms" 1 ), seen 
partly in CI and partly C2, is oriented in a north-south di- 
rection across main disk. Just to the west of this, we identify 
a component blueshifted by 100-200 kms" 1 compared to CI 
that we have assigned to C2. It is not clear whether this is 
a signature of outflow or not. Evidence of split line compo- 
nents in the south-western companion are also seen, with 
velocity separations of 400-500 kms" 1 . This is consistent 
with the fact that ongoing intense star formation is occur- 
ring here, and may imply a superbubble or wind is being 
blown. (Rupke et al. 2005a) find evidence for a neutral gas 
outflow from this system via the Nai absorption line, with 
speeds of ~200 kms" 1 . 



IRAS 20551-4250 

IRAS 20551-4250 comprises of a single nucleus, a relatively 
high surface brightness large-scale tidal tail to the south, 
and a fainter tidal plume to the north. The quasi-elliptical 
light profile of the central regions suggests that the violent 
relaxation of the galaxy collision is almost complete and 
this system is therefore in an advanced coalescence stage. 
Ionized gas emission is only detected in the central regions, 
except for some extended emission to the south-east, at the 
base of the tidal tail. The disk appe ars to be in well-ordered 
rotation with vsini ~ 90 kms" 1 . Mihos fe Bothunj ( |l998| ) 
suggest that these low rotation velocities indicate that we 
are viewing the disk largely face-on. All three line ratios are 
depressed in the central region, as we see in many other 
systems in our sample. 

Much evidence points to this system being at the end 
of the ULIRG phase: a regular velocity field in the disk, cen- 
trally concentrated Ha, and the fact that population synthe- 
sis modelling ( Johansson|1991| suggests that the peak in star 
formation already occurred some 10 7 yr ago. Point-like hard 
X-ray sources have been detected with both XMM-Newton 
dFranceschini et al.|[20Q3| ) and Chandra |Ptak et al.||2003| ), 
and those, together with more analyses of Spitzer /IRS spec- 
troscopy suggest that this system hosts a highly obscured 
AGN | |Farrah et al.||2007| [Sani fc NardmT||2012| , again con- 
sistent with the evolved nature of this ULIRG. In our Ha 
maps, we find evidence for fast-flowing gas that is highly 
suggestive of a turbulent outflow. In a strip across the north- 
ern side of the disk, Ha line widths range between FWHM 
= 500-800 kms" 1 , and the region of peak line widths is 
also blueshifted by 150 kms" 1 beyond that expected from 
pure rotation. The fact that the starburst is no longer ac- 
tive, and the presence of the AGN, strongly implies that the 
outflow is AGN-powered (despite the velocities being quite 
low). Indeed, a fast molecular-phase outflow is identified in 
Herschel/P ACS spectroscopy (E. Sturm, priv. comm), sup- 
porting the AGN-driven scenario. 

The total exposure time allocated for this system was 
split into two, with a second pointing covering the southern 
tail obtained with an exposure time of 6000 s. However, no 
emission lines were detected in the tail region, even after 
binning 5x5 in the spatial axes. 
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IRAS 21504-0628 

This object has received very little previous attention in 
the literature and has never been observed with HST. It is 
classified as a Sy2 although this is uncertain \ Far rah et al.| 
|2QQ3). In optical continuum light, IRAS 21504-0628 appears 
to have a flattened disk with a very faint tidal arm extending 
to the north-west and a small extension to the south-west 



(Due et al. 1997). Our Ha; velocity map indicates ordered 
disk rotation with vsini ~ 175 kms -1 , although there are 
some obvious distortions in the velocity field. The Ha emis- 
sion in the central disk is much less extended, appearing 
compact and symmetrical, although both the northern and 
southern tidal features are also identifiable. The [Sii]/Ha 
and [Ol]/Ha line ratios are depressed in the nuclear region, 
whereas the [Nil] /Ha ratio remains fairly high across the 
whole disk. This is similar behaviour to other galaxies in our 
sample. Gas densities are also measurably high in the cen- 
tral disk. Ha line widths in CI are generally broad (FWHM 
= 250-350 km s _1 ), but in one region to the north of the nu- 
cleus we find line widths up to ^500 kms -1 where we also 
see elevated line ratios in all three indicators. In the central 
region, we identify a cluster of spaxels in which a second (an 
in the very centre a third) Ha component can be fitted. This 
second component is consistently blueshifted by ^50 kms -1 
with respect to CI, but with line widths of <100 kms -1 and 
average line ratios. Although the kinematics and gas proper- 
ties in this system are clearly complex, none of our findings 
is obviously suggestive of an outflow. 



IRAS 22491-1808 

IRAS 22491-1808 has a complex morphology with tidal fea- 
tures extending towards the north-west, the east and the 
south. Two nuclei separated by 2" (~3 kpc) were distin- 



guished in ground-based K-band imaging ( Carico et al. 1990 
Surace et al.|2000| ). However, [Cui et al.| ( [2QQ1| ) find three pu- 



tative nuclei from HST I-band imaging, one located 1.2 kpc 
to the southeast of the central bright nucleus and one 2.4 kpc 
to the west, and therefore suggest a multiple merger sce- 
nario for this system. The many star forming knots found 
around the central region (preferentially distributed around 
the western nucleus) contribute significantly to the total 
bolometric luminosity of the galaxy ( |Surace et al.|20 00). The 
eastern nucleus has very warm optical/near-IR colours, sug- 
gesting that this nucleus may harbour an AGN ( Surace et al. 
2000) . Indeed, [Farrah et al.| ( |2003| find evidence in their SED 
fitting for both a SB and AGN component, where the AGN 
contributes >50% of total IR luminosity. Francesc hini et aT] 
(2003), however, find only very weak X-ray emission from 
this system. Star forming knots are also clearly embedded 
in the tails. 

Our Ha maps show distorted velocities within the nu- 
clear region, as expected from such a complex merger. The 
north-western tail is Ha-bright and shows smooth redward 
velocity gradient. The southern extension is visible in our 
Ha intensity map, although the eastern tail is not. Broad 
C2 line widths (FWHM = 300-350 kms -1 ) are found at the 
base of the north-western tail, suggesting that material here 
has been compressed during the collision and is now form- 
ing stars. The presence of multiple line components in the 
central regions reflects the complex dynamical state of this 



system. As we have seen in many of the other objects in our 
study, the [S n]/Ha and [O i]/Ha line ratios are depressed in 
the disk but rise on the outskirts. The [Nil] /Ha ratio does 
not follow this trend, and is high on eastern side and low in 
north-western tail. 



IRAS 23128-5919 

This system consists of two merging galaxies, the nuclei 
of w hich are separated by a projected distance of 4 kpc 

[Due et al. 



(5"; Zenner & Lenzen 



1993 



1997 



Charman- 



daris et al.|2002 ). As a result of the collision, two tidal tails 



stretch 40 kpc in opposite directions (Bergvall & Johans- 
son|1 985; Mihos & Bothun 1998 ). Based on optical studies, 
the northern galaxy is classified as a starburst, while it is 
unclear whether the southern one is a Seyfert, a starburst 
or a LINER (|Duc et al.||1997l. In hard X-rays the south- 



ern nucleus is detected as a relatively luminous point-like 
source therefore indicating the presence of an obscured AGN 
dFranceschini et al.|2003| |Ptak et aT1|2003| ). 

Previous optical spectroscopy of the southern nucleus 
has shown asymmetric emission line profiles with broad blue 
wings in a region extending ^5 kpc out from the nucleus 
(Bergval l &; Johansson |[T985] ) . Ha emission associated with 
the southern galaxy is significantly extended in the east- 
west direction compared to the I-band morphology. Since 
spectral features characteristic of Wolf-Rayet stars are also 
present over the whole central region ( Johansson & Bergvall 



1988), these observations have been interpreted as evidence 



for a starburst driven outflow from the southern nucleus. 



In their Fabry-Perot observations, Mihos & Bothun (1998) 



also found complex Ha emission line profiles in the southern 
galaxy of this pair, with asymmetric shapes and double- 
line profiles. However, they interpreted the large east- west 
velocity gradient associated with the southern nucleus as 
resulting from a disk seen nearly edge-on. 

These previous studies, however, did not decompose 
the emission line profiles into individual components. In our 
data, the east-west extension associated with the southern 
nucleus seen in the Ha intensity map is characterised by a 
very broad (FWHM = 600-700 kms -1 ) Ha emission com- 
ponent (in our maps assigned to CI). Like Mihos & Bothun 
(1998), we also find this component exhibits a large east- 
west velocity gradient spanning -900 kms -1 . The corre- 
sponding [Sn]/Ha and [Oij/Ha line ratios in this strip are 
low compared to the surrounding material. In the region 
surrounding the southern nucleus, we also see two narrow 
components (in our maps assigned to C2 and C3). C2 is at 
the systemic velocity of the system, while C3 (and CI) are 
blueshifted. The region in which we see three line compo- 
nents is indicated on the Ha C2 radial velocity map with a 
black ellipse. 

Is the region with the large velocity gradient a bipolar 
wind or a disk? It cannot be a significant disk of stars as 
it is not visible on the HST F814W image, or on older op- 
tical continuum imaging (Bergvall & Johansson 1985; Due 
et al.||1997 ). It is also very improbable that a large, fast ro- 
tating, disk of ionized gas is present, without at least some 
associated stars visible in the continuum image. Therefore, 
we conclude that this is a large bipolar outflow. Evidence in 
support of this is: (1) large emission line widths coincident 
with (2) a 900 kms -1 velocity gradient; (3) a lack of contin- 
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uum emission indicating a stellar disk; and (4) the presence 
of a narrow line component at the systemic velocity cover- 
ing the whole inner part of the large velocity gradient region 
(assigned in our maps to C2) that presumably originates in 
the diffuse ionized gas associated with the main body of the 
southern system. Under the outflow interpretation, we are 
seeing both the red and blueshifted plume of a large-scale 
bipolar wind with outflow velocities of >450 kins" 1 (given 
the unknown inclination). The Av max of 800-1000 kms -1 
(Table [3]) and the known presence of an AGN in the south- 
ern nucleus together suggest the outflow is AGN driven. In- 
terestingly, no evidence for a molecular gas outflow is seen 
in Herschel OH line spectroscopy (E. Sturm, priv. comm.). 

The aforementioned narrow line component seen in C2 
at the location of the outflow plumes matches consistently 
(in terms of line width and radial velocity) with the CI com- 
ponent seen away from the bipolar wind. We therefore asso- 
ciate the two as originating in the same material (and were 
just assigned to the different components as a result of our 
component assignment convention; see Section 2.2). Under 
this assumption, the ionized gas not associated with outflow 
plumes shows no evidence of rotation (even after reexamin- 
ing the velocity map with a —100 — hl00 km s _1 scaling). The 
[Nil] /Ha and [Sii]/Ha line ratios for this systemic velocity 
component are high, consistent with shock or AGN excita- 
tion, and a depression is seen in the [Sll]/Ha and [Oij/Ha 
ratios at the location of both nuclei. Both nuclei also show 
up in the electron density map with peaks at 400-500 cm -3 . 

The proximity of the nuclei and complex dynamics led 



Mihos & Bothun (1998) to suggest "we are most likely wit- 



nessing IRAS 23128-5919 in a very short-lived phase where 
violent relaxation is rapidly altering the structural and dy- 
namical properties of the system" . This system is therefore 
very worthy of detailed follow-up observations. 



